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Summary: An attempt is made to assess the rate at which failures occur within a 
group of similar structures, or structural components, when each structure is subjected 
to the same set of applied loads. The structures are similar in every respect except for 
small inherent differences in strength from one structure to another. On infrequent 
occasions the applied loads may be so intense as to destroy the weaker members of 
the group. For the problem in which the structural strengths are distributed within 
he group in a “ normal” fashion, and when the frequency of severe loads diminishes 
exponentially, it is found that the rate at which failures occur within the group of 
structures is governed by two simple non-dimensional parameters. Numerical values 

of these parameters are studied more fully. 


1. Introduction 


In the design of engineering structures it is assumed generally that all com- 
ponents of the same type, and manufactured from the same material, have equal 
strengths. It is well known, however, that all constructional materials—whether 
produced naturally or manufactured—show distinct variations in strength from one 
est sample to another, even when special efforts are made to produce a material of 
onsistent strength. With a set of structural components manufactured from a 
atural material, such as timber, the variations in strength from one component to 
unother may be considerable; those materials which are manufactured—such as steel, 
bluminium alloy and, to a limited extent, concrete—may be controlled in quality 
more carefully, and for these materials it is possible to prescribe certain minimum 
strengths with reasonable accuracy. However, even for the manufactured materials 
e testing of a number of structural components reveals inevitable variations in 
strength from one component to another, although the range of strength scatter 
ay be relatively small. In addition to inherent material variations, there is the 
urther difficulty that manufacturing processes cannot provide structural components 
which are exactly similar in geometrical form, in surface finish, and so on. It may 
concluded that, although the mean strength of a number of similar structural 
omponents is usually defined closely, the scatter of strength among those compon- 
nts is governed by the quality of the material and the extent to which the 
anufacture of the components is controlled. 


A further important consideration in assessing the safety of a group of similar 
uctures is that extremely severe loads—more intense than the design loads—may 

be encountered by a few of the structures during their working lives. The infrequent 
urrence and uncertain nature of these loads may tempt the designer to ignore 
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their effects, because the life expected of each structure may be so short that there 
is little chance of severe loads being encountered. For example, hurricanes and 
earthquakes may destroy completely certain light types of building; obviously it is 
not economical to design the less important structures to withstand the most severe 
effects of such disturbances. To take account of these severe loads and to provide 
the structures with the necessary additional strength, would mean that a longer life 
is expected of each structure, and in many instances this is unjustifiable. A more 
logical approach to this problem is to accept an occasional failure of one of a group 
of structures, and at the same time to provide the structures with sufficient strength 
to give a frequency of failures which is not unreasonable for structures of that 
particular type. In the aeronautical field the design of aircraft wings involves a 
similar problem; on rare occasions an aircraft may encounter extremely severe 
upward or downward gusts. The frequency of occurrence of gusts of such intensities 
is so small that it is not economical to encumber every aircraft with the additional 
structural weight necessary to provide against these severe loads. 


On considering these two aspects of the problem, only qualitative conclusions 
can be deduced about their relative importance at this stage. However, a discussion 
of the frequency of failures among a group of similar structures will revolve 
obviously around considerations of 


(i) the variation of strength within the group of structures, and 


(ii) the frequency of occurrence of severe loads. 


If the strength variation is considerable, then it may be found that both moderate 
and severe loads destroy some of the structures. On the other hand, if the strength 
variation is small, then only the severest loads cause destruction among the group 
of structures. 


The first writer to study the problem in simple terms was Pugsley“’ who 
discussed the case of a group of aircraft wings of variable strength, each wing 
encountering gust loads of varying intensities. The approach used by Pugsley is 
extended in the present paper to give finally a method of analysis which is wider 
in scope. 


Notation 
number of structures within the group 
intensity of load applied to a structure 
number of structures of strength less than P 


number of loads of intensity greater than P applied to each structure 
in unit time 


number of failures in N units of time 

frequency of failures 

mean strength of the group 

standard deviation of strengths within the group 
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INCREMENT IN P 


Fel R 


LOAD INTERVALS 
Fig. 1. 
Distribution of strength and variation of applied loads on the basis of finite differences. 


y normal distribution function 
k aconstant related to the applied loads 
(x) standardised normal distribution function 
(x) standardised normal cumulative distribution function 
non-dimensional strength parameter (=o /P,,) 
8 non-dimensional load parameter (= KkP,,) 


2. Solution in Terms of Finite Differences 


2.1. Distribution of strength and variation of applied loads 

Consider a group of N structures, or structural components, similar in all 
respects except that the strengths may vary from one structure to another. Suppose 
A, structures have strengths less than that corresponding to P, (Fig. 1); for con- 
venience the load scale is divided into a number of equal finite intervals. If the 
structure is a simple tension member of a framework, then all structures would be 
expected to withstand loads greater than zero load; in this case, A,=0. Also there 
must be an upper limit for the strength, and it is assumed that all structures are 
weaker than some finite strength Pz; then A,=N. 
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The number of loads applied to each structure in unit time, and in excess of 
P,, is taken as E,. It is assumed that each structure of the group is called upon to 
withstand the same set of loads. The unit of time is defined quite arbitrarily; in 
some problems, such as those associated with aircraft wings, the unit of time is 
replaced for convenience by a suitable unit of distance. The time unit is simply the 
most convenient basis for comparing the relative frequencies of the applied loads. 
Since all the applied loads must be greater than zero intensity, E, represents the total 
number of loads applied to a structure in unit time. The variation of the function 
E, must be derived from the measurement of loads in actual structures of the type 
under consideration. 


2.2. Increments in strength (Pugsley’s method) 


Consider the number of structures of strength P,; there are A, structures of 
Strength less than P,, and A,,, structures of strength less than P,,,. So that 
(A,4,—A,) structures have strengths within the interval P,; these structures are 
destroyed only by loads which exceed P, in intensity, that is, by E, loads in unit 
time for each structure. If it is assumed* that whenever structures of a certain 
strength are destroyed then they are replaced by new structures of the same strength, 
the effects of the loads on the structures can be imagined as a continuous process of 
destruction and replenishment at constant strength. The number of failures suffered 
by structures of strength P, is E,(A,,,—A,); so that the total number of failures 
suffered by structures of all strengths is 

r=R-1 
since A,=0, and (A,,,—A,)=0 for r= R; the failures F, occur in N units of time 
Equation (1) is a formal statement of the method used by Pugsley". 


2.3. Increments in applied load (alternative method) 

There is an alternative method of deriving the total number of failures suffered 
by the structures; consider, for example, the damage caused by the loads of intensity 
P,. The number of loads applied to each structure in unit time and in excess of P,_, 
is E,_,, while the number of loads in excess of P, is E,. So that in the interval P, 
there are (E,_,— E,) applied loads, which destroy structures of strength less than P.. 
On the basis of replenishment at constant strength, the number of failures caused by 
these loads is , (E,_,— E,). The total number of failures caused by all applied loads 
is therefore 


1 
A, —E,)+ArEr_,. (2) 


=¢ r=R 


r 


1 


since A,=0, Ap=Ans, =Any.=----=A,, and E,, is assumed equal to zero. The 
failures F, occur again in N units of time. 


*This assumption is considered by Pugsley") to be a reasonable approximation to practical 
conditions. 
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2.4. Comparison of the two methods 

Both methods give values of the total number of failures occurring within the 
group of N structures. The former method (Pugsley’s approach) gives the number 
of failures suffered by the group of structures, while the latter method gives the 
number of failures caused by the applied loads. Since the total number of failures 
in each case is the same, F, and F, are equal; to demonstrate this, equation (2) is 
written first in its expanded form 


F,=A, (E,-—E,)+A, (E, ~E,)+. En-,) + AnEn-). 
The terms in this series may be rearranged to give 


F,=),E,+E, (A,—A,)+E, (A, —A,)+ Ar- Xr-,), 


r=R-1 


so that F,=E,\,+ 
r=1 


TABLE I 


ESTIMATION OF THE FREQUENCY OF FAILURES AMONG A GROUP OF STRUCTURES 
A,=number of structures of strength less than P,, 
E,=number of loads in excess of P, applied to each structure in unit time. 


Increments in X Increments in E 
Ar EAX AE AAE 


Interval 


NRO 


Total number of failures=60,385 in 10* time units 
=60+ in 1 time unit. 
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1 & 15 400 40 200 8,000 23 9,200 
P= 16 600 25 176 4,400 15 9,000 
p 17 776 16 121 1,936 6,984 
i. : 18 897 10 65 650 5,382 
y § 19 962 6 27 162 3,848 
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24 1,000 0 0 0 1,000 
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14,000 
[| NUMBER OF POSITIONS OF PEAKS CALCULATED 
FAILURES 4" FROM EQUATIONS (27) AND (28) 
12,000} / 
10,000; INCREMENTS 
INA 
8,000} 
6,000} 
4,000} 
2,000}; 
LOAD RANGES 
Fig. 2. 


Distribution of failures for the example of Table I. 


Although the summations of the two series are the same, corresponding terms 
in each series are not equal. The sum of terms of equation (1) up to r=q is 


T=q 
F, (q)=E,A, = E, —A,). 


From equation (2), the sum of terms up to r=q is 


Tr 


=4q 
F, (q)= Xr, ~E,), [q < R]. 


=] 
so that F, (q)—F.(q)=E (3) 


It is deduced that F, (q)> F, (q); at the low end of the strength range and the upper 
end of the applied load range the two functions are equal. 


The two methods are compared in the solution of the frequency of failures for 
the group of structures outlined in Table I. The group consists of 1,000 structures 
with a symmetrical distribution of strengths about a mean value of 15. The variation 
of applied loads is not necessarily typical of any practical problem; it is used in this 
example to illustrate the important differences between the two methods of solution. 
The individual terms of the two series given by equations (1) and (2) are plotted in 
Fig. 2. The largest term in equation (1) corresponds to those structures suffering the 
greatest damage; in this light Pugsley’s method is useful in indicating the weakest 
range of strength. On the other hand the largest term of equation (2) corresponds 
to those applied loads causing the greatest failure frequency; equation (2) indicates. 
therefore, the range of loads causing the greatest damage. 


The examples discussed by Pugsley“’ may be analysed by either equation (1) or 
equation (2); although the total summations given by these equations are the same. 
there is again a marked difference between the positions of the largest terms in the 
two series, 
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3. The Use of the Differential Calculus 


The main purpose of the present paper is to determine those factors which 
govern the frequency of failures; it is advantageous to simplify the preceding analysis 
by replacing the finite intervals of load and strength by infinitesimal increments. 
This will give continuous functions for strength distribution and applied 
load variation. 


Suppose E(P) represents the number of times the applied load P is exceeded 
in intensity on each structure in unit time; suppose also that in the group of N 
structures there are A(P) having strengths less than P. The functions E and X have 
the continuous forms shown in Fig. 3; E, defines again the total number of loads 
applied to each structure in unit time; it is assumed that E—>0 as P—>oo. The 
total variation of strengths is finite, so that 4(0)=0 and A(oo)=N; that is, all 
structures are stronger than zero strength and weaker than infinite strength. 


ns 


Fig. 3. 

3) Continuous variations in E and A, 
er 

3.1. Infinitesimal increments in strength 
or Of the N structures in the group, 6A have strengths in the range from P to 
as (P+é6P). Each of these structures withstands E loads in excess of P in unit time. 
yn On assuming as before that a destroyed structure is replaced by a structure of the 
is same strength, the number of failures suffered by structures in the strength range 
n. 5P is Edd. The total number of failures suffered by all structures is 
in 
st = | [ (in N time units), . ‘ (4) 
1s P=0 0 
S, 

where A’=dA/dP. The number of failures in unit time, or the frequency of failures, 

is therefore 
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3.2. Infinitesimal increments in applied load 


Again, consider the applied loads in the range from P to (P+ 45P); there are E 
loads in excess of P, and (E+ 5E) loads in excess of (P+ 45P). So that in the range 5P 
there are —5E applied loads. The number of structures destroyed by the —45E loads 
is —AdE, since A structures are weaker than P. On the basis of replenishment at 
constant strength the total destruction due to all applied loads is 


| —\A\dE=- | En’dP, (in N time units), ; (6) 


0 
where E’=dE/dP. The number of failures in each unit of time is 


1 


3.3. Comparison of the two frequencies of failure 


On integrating equation (6) by parts, 


But E(0o)—0, and (0)=0; so that [ae | 


F,= ENdP=F,. 
0 


As for the finite difference method, the total number of failures suffered by all 
structures may be derived by either equation (5) or equation (7). 


On the basis of increments in strength, the total number of failures suffered by 
all structures of strength less than P is 


F,(P)= { evap. ‘ (8) 


Alternatively, for increments in applied load, the total number of failures caused by 
loads less than P is 


F, (P)= - AE’AP. 
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Equation (9) may be written 


P 
0 


F,(P)=— { evar] 
since A(0)=0. 


It is concluded that F, (P) is always greater than F,, (P), unless A or E is zero. 


4. Idealised Variations of Strength and Applied Load 


4.1. Strength variations 


The analysis of Section 3 can be used most profitably if it is possible to define 
more precisely the forms of strength distribution on the one hand, and the variation 
of applied loads on the other. 


For the distribution of strengths of the structures a normal distribution about 
the mean strength is assumed. This assumption disregards any deviations from 
symmetry which actual strength variations may show. On this basis, the number of 
structures having strengths in the range from P to (P+4P) is y4P, where y is the 
normal distribution function defined by 


N (P— Pn)? 


in which N=total number of structures in the group, 
P,,=mean strength of the group, and 
o=standard deviation of the strengths of the group. 


A(P) has been defined as the number of structures of strength less than P; so that 


Then A is the normal cumulative distribution function. 


4.2. Load variations 

For the definition of the applied loads it is assumed that severe loads occur far 
less frequently than more moderate loads; extensive studies of gust loads on 
aircraft wings suggest that the function E is defined closely by the relation 


in which E(P)=total number of loads in excess of P in unit time, 
E,=total number of loads in excess of zero load, and 

k=a positive constant [having the dimensions of (load)~"]. 
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Equation (13) is by no means true of all structural loading conditions; for example, it 
has been suggested’ that wind loads on buildings are distributed normally about a 
mean wind velocity of approximately 23 m.p.h.; problems such as this require 
special study. 


4.3. The frequency of failures 

The idealised strength distribution and applied load variations are shown in 
Fig. 4; the strength function A is anti-symmetrical about the mean strength P,,. 
Hypothetically the total range of strength scatter is from P= —0o to P=+0. 


> 
Fig. 4. 
Idealised representation of the strength distribution and applied load variations. 
The variation of positive loads is given by equation (13); suppose, for con- 
venience, that all loads are negative initially, and increase progressively up to their 
maximum values; then for P< 0, we have E(P)=E,,. 


In terms of increments in the strength function A, the total number of 
failures is 


+% 0 x 
F,= [ EN dP=E, | vaP+ | Exar. (in N time units), 


The frequency of failures is 


f= + EN | , (in unit time),. . (14) 


where A, is the value of A at P=0O. 


On the basis of increments in the applied load function E, the frequency of 
failures is 


+X 


~ AE’dP= N 


ox 


AE’dP, (in unit time), . 


j 
0 


since E’ (P)=0 for P < 0. 


E 
7 
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Equation (15) may be written 


On substituting for A and E the values given by equations (11) and (13), 


= lave exp { - } dP, E=E, exp(— kP). 


Then 
f= | exp { - } dP + J expc- kP) exp { aP | 


which may be written 
exp(- x?/2) dx-+exp (47k? kP,) | exp(—x"/2)dx 
(ox- 


(18) 


Now, — exp(- x?/2)dx is the standardised normal cumulative distribution 


function, written here for convenience as (x). Then equation (18) takes the form 


- fs) + exp (Jo?k? — a ~ok), (19) 


since exp(—x?/2)dx= exp (—x*/2)dx= (x). 


When the first term, ®(—P,,/o), in equation (19) becomes large in relation to the 
second term, the idealised representation of the frequency of failures loses all 
practical significance, because the weaker members of the group, for which P <0, 
become important. A study of the relative values of the two terms in equation (19) 
indicates the stage at which the idealised problem is meaningless. 


4.4. Distribution of failures 
The failures suffered by structures up to strength P is 


P 


F, (P)= EN AP. 
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We are interested in values of P > 0; so that 
0 P P 
F.@M=E, + ( EN’ dP=E,A,+ | EN 
0 0 


Then dF . . .  . (20) 


where y is the normal distribution function. On substituting for y, 


F,’(P)= exp (407k? — kP,,) exp { (Pa } (21) 


so that F,’(P) is a type of normal distribution function with values disposed sym- 
metrically about a mean value P=P,,—o°k; the maximum value of F,’ (P) is 


Phun ™ exp (40°k? — kP,,,). 


It may be concluded that the most heavily damaged strength range occurs at 
P=P,,—o°k, that is, at a value of P below the mean strength P,, of the structures. 


The total number of failures caused by loads up to an intensity P is 


F,(P)= - | aap. 
We are interested in values of P>0; so that 
> 
F,(P)=- 


and dF (P)/dP=F,’ (P)= —AE’. 


On substituting for A and E 


a 


F,’ (P)=E,k exp (— kP) [ + 


2c? 


The maximum value of F.’(P) occurs when F.” (P)=0. Now 


F,” (P)= + E,Nk exp (~ kP) oko (F Pn) | 


where #(x)= the standardised normal distribution function. 
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Then F.,” (P) is zero when 


PP, . ‘ ‘ (24) 


The required value of P is the root of equation (24), which cannot be solved exactly 
in a simple fashion. It is interesting to note that F,,’ (P), unlike the function F,,’ (P), 
is not a normal distribution. 


5. Basic Parameters of the Idealised Problem 


The important feature which emerges from the idealised representation of the 
strength distribution and the applied load variation is that the frequency of failures, 
as given by equation (19), is governed by the two parameters (c/P,,) and (KkP,,). 


Suppose kP.. = B. ‘ 


Then the frequency of failures is given by 


The function (f, /E,) defines the number of failures in unit time for each load applied 
to a structure; 2 is a non-dimensional parameter which gives a measure of the 
scatter of strength within the group of structures; the parameter 8 defines the 
frequency with which the applied load P,,, is exceeded. If there is no scatter of 
strength then 1/2=00, and 


(f,/E,)=® (— 00) + exp (— 8) (+ 00)=exp (— A). 


In this case, f, =E, exp(— 8)=E, exp (— kP,,,), which is the frequency with which P,, 
is exceeded; of course, this would be expected to correspond to the frequency of 
failures when all the structures have the same strength P,,,. 


In terms of « and 8, the most heavily damaged strength range is given by 


Moreover, the most heavily damaging load range is given by the solution of the 
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ort 


005 O15 020 025 
VALUES OF 


Fig. 5. 
Chart of idealised frequencies of failure. Points X,, X,, X, refer to Pugsley). 


Since the function f,/E, is governed by only two parameters according to 
equation (26), values of f,/E, may be studied on a single chart, as indicated in 
Fig. 5. In this diagram the failure frequencies are given for values of the “strength” 
parameter «/P,, from 0:00 to 0:25, and for values of the “applied load” parameter 
kP,, from 0 to 50. Low values of kP,,, have no practical significance, but are included 
for completeness. In many problems o/P,, is very small; Fig. 5 shows clearly the 
extent to which the failure rate is increased by the presence of even small scatters 
of strength; for example, if AP,,—40, an increase of the strength parameter o/P,, 
from 0 to 0-05 produces a tenfold increase in the failure rate. 


The three numerical problems discussed by Pugsley? are indicated in Fig. 5 by 
the points X,, X, and X,; in the first example the strength properties of the structures 
give «/P,,=0-097, while the applied loading conditions approximate to kP,, = 15:0. 
These values of the two parameters correspond to the point X, in Fig. 5; the fre- 
quency of failures given by Fig. 5 is 31 per unit time, compared with 32 per unit 
time computed from Pugsley’s finite difference formula. In the second example the 
mean strength of the group of structures is reduced to give 7 /P,, =0-115; the value of 
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oF (ir) kP,=0 
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025 0-05 0-10 


VALUES OF (/r,,) 


Fig. 6. 
Position of the most heavily damaged strength range. 


kP,, is reduced to 12-7; these correspond to the point X, of Fig. 5. There is 
obviously a concomitant increase in the frequency of failures. In the third example 
Pugsley considers the effect of reducing the standard deviation of the strength distri- 
bution; in this case /P,,=0-049 and kP,,= 15-0, giving the point X, in Fig. 5. The 
change in failure frequency is not so pronounced because the change from X, 
to X, takes place along the flat portion of the kP,, contour. 


From equation (27) can be derived the most heavily damaged strength range in 
terms of 7/P,, and kP,,; Fig. 6 gives values of P/P,, which satisfy equation (27). 
In all practical problems the most heavily damaged strength range is below the mean 
strength P,,; this can be checked for the three examples discussed by Pugsley“’. For 
the most heavily damaging applied loads equation (28) is solved for P/P,,; the 
results of a number of computations are given in Fig. 7. For the numerical example 
studied in Table I, the positions of the turning points of the two curves, as 
calculated from equations (27) and (28), are indicated in Fig. 2; in that example the 
variation of applied loads follows the exponential law of equation (13). The 
estimated positions of the maxima are in excellent agreement with the actual 
turning points. 


6. Conclusions 


The solution of the problem in terms of infinitesimal increments in strength and 
applied load leads to results which may be used to considerable advantage in study- 
ing the problem more generally. The frequency of failures can be derived in two 
ways—the one based on strength increments and the other on applied load 
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VALUES oF (7/p, ) 


0°4 


0-2 


0-05 


Fig. 7. 
Position of the most heavily damaging range of applied loads. 


increments; the result is the same in each case, although there are different 
distributions of failures among the load ranges. The one method indicates the most 
heavily damaged range of strengths, and hence the weakest section within the group 
of structures, while the other method shows clearly the most damaging range of 
applied loads. 


For the idealised variations of strength and applied load, the frequency of 
failures is governed by “strength” and “applied load” parameters; in this respect 
Fig. 5 is useful in indicating changes in the frequency of failures which accompany 
alterations of either the mean strength or the scatter of strength of the group 
of structures. 
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The Effect of Delta Vanes on Supersonic Wings 


M. A. GORGUI 
(Department of Mathematics, University of Manchester*) 


Summary: For supersonic wings it is known that the presence of a tip reduces the 
lift and the performance compared with the corresponding values for the two- 
dimensional wing. It appeared from some work by K. Stewartson® that the addition 
of a vane at the wing tip might counteract this effect. The purpose of this paper is to 
determine the effect due to the presence of the vane. For this purpose a wing in the 
form of a semi-infinite wedge whose leading edge is swept back is supposed placed at 
incidence to the flow, with a symmetric triangular vane with diamond-shaped base 
attached to it, its vertex being at the tip of the leading edge, and its plane of symmetry 
parallel to the flow. The problem is of the “ conical field” type and is solved by the 
method of S. Goldstein and G. N. Ward. 


The problem is divided into its anti-symmetrical and symmetrical parts. The 
results obtained in Ref. 1 are applied to determine the field in each case, both outside 
and inside the Mach cone of the vertex. 


In the anti-symmetrical part, three cases are considered, with a different vane in 
each case, and the lift coefficient for the parts of the wing inside the Mach conc is 
determined. In the symmetrical part, only the case when the vane lies entirely inside 
the Mach cone is considered, and the drag force is obtained. The performance is then 
calculated, and curves are drawn to show how it is affected by the existence of the vane. 


1. Introduction 


In supersonic flow the lifting force at the tip of a finite plane wing with a 
straight leading edge is less than the value for a two-dimensiona! wing with the same 
angle of sweep of the leading edge. The two-dimensional value of the lift is attained 
at points on the wing which lie just outside the Mach cone with vertex at the tip of 
the leading edge. If the tip is straight and parallel to the undisturbed stream, then 
the lift falls from the two-dimensional value at the Mach cone to zero at the tip. 
The drag changes in a similar way but much less rapidly; it decreases from the 
corresponding two-dimensional value at the Mach cone to a value that depends on 
the angle of sweep. It seems then that the performance of such a wing is not as good 
as that of the corresponding two-dimensional wing. But from some results obtained 
by K. Stewartson it appeared that if a vane is attached to the wing tip, it might 
counteract this effect. For the reduction in lift is due to flow around the tip from 
the pressure to the suction side, and this would presumably be reduced by the 
presence of a vane. These considerations suggested the problem of studying how 
the addition of a vane to a wing will affect the lift, the drag and hence 
the performance. 


For this purpose, the wing is taken in the form of a semi-infinite wedge, whose 
Straight leading edge is swept back and whose tip is along the undisturbed stream. 


*Now Department of Mathematics, Alexandria University. 
Originally received December 1952. 
[The Aeronautical Quarterly, Vol. V, November 1954] 
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Fig. 1. 


The vane has a triangular plan form and a diamond-shaped base; it is attached to 
the wing tip such that its vertex is at the tip of the leading edge, and their planes of 
symmetry are perpendicular to each other. The whole is placed in the flow with the 
wing at incidence 2, and the vane at zero incidence. A sketch of the wing and 
the vane when it lies entirely inside the Mach cone is shown in Fig. 1. 


The edges of the wing and the vane are taken straight, so that the problem is of 
the conical flow type, and it will be assumed that the incidence of the wing and the 
thickness of both the wing and the vane are small, so that the linearised theory of 
conical fields as set out in Ref. 1 can be used. Moreover, for the purpose of calcula- 
tion, the wing and the vane are supposed to extend to infinity downstream. If the 
wing and the vane have supersonic trailing edges, the solution thus obtained will be 
valid upstream from the trailing edges, and the solution for a system of finite extent 
in the stream direction is determined. 


The problem is solved by dividing the velocity field into its anti-symmetrical 
and symmetrical parts and treating each independently; the results may then be 
superposed to give the complete solution. The anti-symmetrical part corresponds 
to a wing of infinitesimal thickness at incidence 2 with a plane vane normal to the 
wing, and from this solution the lift is calculated. The symmetrical part corresponds 
to a wing and vane of small finite thickness, of the shape described earlier, but at 
zero incidence. This solution will give no lift, but only a drag force. 


The results of the calculation show that the effect of the vane on the lift is to 
counteract the loss of lift caused by the presence of the wing tip. The magnitude of 
this effect depends on the angle of sweepback as well as the vertex angle of the vane. 
This effect is shown in Figs. 7 and 8, and the formula for the lift coefficient when 
the trailing edges are perpendicular to the stream is given by equations (60), (61) and 
(62). The formula for the drag is obtained in equation (93), and its variation with 
the vertex angle of the vane is shown in Fig. 10. 


An index of merit for the performance (the lift/drag ratio) is then calculated, 
and it is found that it first increases and then decreases as the semi-angle of the 
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vertex of the vane increases up to the Mach angle. The maximum performance for 
a particular vane is found to increase with the decrease of the vane thickness. It 
appears that a wing-tip vane, if its dimensions are suitably chosen, could be used 
in practice to improve the performance of a supersonic wing. 


Notation 
a constant defined in equation (20) 
b constant defined in equation (23) 
2c length of the cross section of the vane by the plane z=cot p 
f(p) arbitrary periodic function of «+6 
k,k’ modulus and complementary modulus of elliptic functions in Section 6 
l,m,n — integers in Section 4 


pa variable of integration in equation (8), and the degree of the function 
H (z) in equations (27) and (65) 


polar co-ordinates in the x,y,-plane 
t complex variable defined in equation (17) 
t’ imaginary part of ¢ 
real values of t 
u,v,w perturbation velocity components 
x,y,Z rectangular Cartesian co-ordinates 
X,,y¥, | rectangular co-ordinates in the plane z=cot » 
z complex variable defined in equation (22) | 


5 ¢in Sections 4-7 
Z value of z defined in equation (23) J 


A,B _ real constants defined in equations (30), (29) 
C,C,,C, real constants defined in equations (76), (69) 
lift coefficient 
D total drag 
D, drag over the wing 
D, drag over the vane 
D, drag due to the suction force 


E(x),E incomplete and complete elliptic integrals of the second kind to 
modulus k 


F(,\),F incomplete and complete elliptic integrals of the third kind to 
modulus k 


F(p) arbitrary periodic function of « — 4 in equation (8) 
G,,G,,G, functions of t 
H,,H,,H, functions of z 
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complete elliptic integrals of the first kind to modulus k, k’ 
lift force 
Mach number 
P perturbation pressure 


W velocity of undisturbed stream 


a angle of incidence of the wing 


angle of sweepback of the leading edge of the wing 
y | semi-vertex angle of the vane 
Y. Value of y at which the lift is a maximum 
y, value of y at which the lift is restored to the two-dimensional value 
26 thickness of the wing 
thickness of the vane 
é,,¢, constants defined in Section 7 
tsi values of <, and ¢, respectively in the plane z=cot p 
6 a variable of integration in equation (73) 
A a constant defined in Section 7 
Mach angle 
p density of the undisturbed flow 
6 variables defined in equation (6) 
uv, value of « for the leading edge of the wing 
«, Value of o for the upper tip of the vane 
x complex variable defined in equation (72) 
yY variable of integration in equation (74) 
» performance 
KR “real part of” 


F “imaginary part of” 


2. The Linearised Theory of Conical Fields 


Let (x, y, z) be a system of rectangular co-ordinates with the origin at the vertex 
of the conical field (the tip of the leading edge in this problem) and the z-axis in the 
direction of the undisturbed stream, M its Mach number, and »(=cosec~'! M) its 
Mach angle. The flow is assumed to be irrotational, so that the perturbation velocity 
is the gradient of a potential function ¢. Then to a linear approximation, the 
rectangular components u, v, w of the perturbation velocity and the potential 
satisfy the differential equation 
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and from the linear approximation to Bernoulli’s equation, 


where P is the perturbation pressure, and p is the density of the undisturbed stream. 


New co-ordinates (x,, y,) and (r, 4) are now defined by 


= 6 — 
Th 
(3) 
y,=rsin 6- 


The Mach cone of the vertex then corresponds to the circle r=1 in the x,y,-plane. 


In conical flows the velocity components u, v, w are functions of (x,, y,) or 
(r, 9) only, since they are homogeneous functions of degree zero in (x, y, z). In terms 
of (r, 4) the equations (1) satisfied by u, v and w become 


(r? +r(2r?—1) =| (u,v, w)= (u,v,w). . ‘ (4) 


The theory of such flows is developed in detail in Ref. 1 and such results as are 
needed for the present problem are quoted from this. 


3. The Anti-symmetrical Part 


The wing and the vane in this case are of infinitesimal thickness. The wing is 
at incidence a while the vane is at zero incidence, so that the boundary conditions of 
zero normal velocity on both the upper and lower surfaces of the wing and on both 
sides of the vane, are respectively 


v= —Wa 
(5) 


and u=0. 


The solution of equation (4) satisfying the boundary conditions (5) varies for the 
different values of r. It is therefore necessary to obtain the solution for the cases 
r>1 and r<1 separately. In other words the two fields (outside and inside the 
Mach cone) must be solved separately. 


The field outside the Mach cone of the origin 


The general theory is dealt with in Ref. 1. The x,y,-plane is transformed to 
the o6-plane, defined by 


xX, = sec 7 cos 8, 
(6) 


y, = seco sin 
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Then equation (4) for r> 1 will be transformed into 


0? 0? 


The part of the x,y,-plane for which r> 1 corresponds, in the «@-plane, to a strip 
between and The axis y,=0 corresponds to (for x, >0) and to 
6=n (for x,<0). The characteristics of equation (7) are the two sets of straight 
lines « + 6=constant and «—6=constant. The first set corresponds in the x,y,-plane 
to those parts of the tangents to r= 1 on which @ increases as the circle is approached, 
and the second set corresponds to those parts of the tangents on which @ decreases 
when the circle is approached. The solution of equation (7) for irrotational flow was 
found in Ref. 1 in the form 


o+6 
u= | f (p) cos pdp + | F (p)cos pdp, 


v= J sin pdp+ F(p)sinpdp, (8) 


o+6 


w= —tany f(p) dp —tan nf F (p) dp, 


where f(p) and F(p) are periodic functions of period 27. They are continuous 
everywhere, except perhaps on crossing certain specified characteristics. 


Let «, (shown in Fig. 2) be related to the angle of sweepback / by the relation 
sec 7, =cot » cot (9) 


and let 2c be the length of the cross section of the vane by the plane z=cot 1, so that 


where 2y is the vertex angle of the vane. Three different cases must then be 
distinguished : — 


(1) c<l, the vane lies entirely inside the Mach cone, 


(2) 1<c<cosec«,; in this case the edges of the vane lie in the regions 1 
and 2 in Fig. 2, 


(3) © >coseco,; and the edges of the vane lie in region 7 in Fig. 2. 


Case 1. c<l 


In this case the field outside the Mach cone is not affected by the vane, but by 
the wing alone, and this will affect only the part of the field bounded by the simple 
wave and the shock wave produced by the leading edge of the wing in the upper 
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Fig. 2. 


and lower half space. The approximation used in the linearised theory reduces these 
waves to surfaces of discontinuity represented by the characteristics k, k’ (Fig. 2). 
Hence the velocity components are“ zero everywhere outside the Mach cone, except 
in regions 1 and 2 in Fig. 2, where they are 


Region Uu v Ww | 
1 —Wacotc, —Wa Wacosec o, tan 
Wacotc, -Wa -—Wacoseco,tann |} 


Case 2. 1<c<cosec 


This case is shown in Figs. 3(a) and 3(b). The corresponding regions in the 
x,y,-plane and the o6-plane are numbered to correspond; the tip of the vane in 
Fig. 3(a) corresponds to the point (*/2, in the o@-plane. 


In the region 1 the velocity components are the same as in the first case and are 
ot+d 


given in equation (11). For u to be zero on the vane, f f (p) cos p dp must decrease 

discontinuously by — Wacot«, as (« + 4) decreases through =—«,, i.e. when region 3 

in the «6-plane is entered along « - #= constant; also | F (p) cos pdp must decrease 


discontinuously by — Wacot«r, as «— 6 decreases through o,, i.e. when region 4 in 


the «#-plane is entered along «~+4=constant. Hence for region 3 


The discontinuous decrease in 7 is 
singe, Lim | , 
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and that in w is —Watanucot«c,seco,. Hence 


u=0, 
w= 


throughout region 3. 


The values in region 4 are found similarly to be 


u=0, 
v= —Wa(l-—cote, tanc,), (15) 
w= W2 tan cosec (1 — COS To sec o,). 


and the velocity components vanish everywhere outside these regions. 
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Case 3. c>coseco, 


In this case the simple wave k produced by the wing tip is reflected on the vane 
(see Figs. 4(a) and 4(b)), thus increasing the number of regions to be considered. 
In the same way as before, the velocity components are found to be 


Region u v w 
| -Wacotc, —Wa Wz tan cosec 7, 
2 0 0 0 
. (16) 
3 0 0 0 
4 —-Wacots, Wz tan cosec 
5 0 —2Wa 2Watanpcoseco, 


and outside these regions the velocity components vanish everywhere. 
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4. The Field Inside the Mach Cone of the Vertex 


As in Ref. 1, the velocity components inside the Mach cone (i.e. inside the 
circle r= 1 in the x,y,-plane) are determined by transforming the interior of r=1 into 
a complex t-plane defined by 


cos6—i(1—r’*)' sin 6 
t r 


(17) 


The circle r=1 is transformed with its interior into the whole of the f-plane, 
with cuts along the portions of the real axis for which |t|=1: the upper half 
0<6<* is represented on the upper half-plane, ./ (t)>0, and the lower 
half —-z< 6<0 on the lower half-plane, .7 (t)<0. The axis x,=0 inside the unit 
circle is transformed into the imaginary axis in the t-plane, while the axis y,=0 
inside the unit circle is transformed without any change into the same part of the real 
axis in the f-plane. 


The velocity components u, v, w of the perturbation velocity are found to 
satisfy Laplace’s equations in this t-plane and therefore 


u= {G, (t)}. | 
v= R{G.(}, (18) 

w= & {G.(}, | 
where G, (1), G.(0), G,(t) are functions of the complex variable ¢; from the 
condition of irrotational flow 


G,’(t): G,’ : —ttanp. » 


To obtain the field inside the Mach cone it is necessary to find the complex 
functions G, (t), G.(t), G,(t), which are related by equation (19) and whose real 
parts satisfy the boundary conditions on the wing, the vane, and the Mach circle. 
Since the boundary conditions on the Mach circle differ for the different cases con- 
sidered before, these functions must be found for each case separately. 


Case 1. c<l 


The projection of the system in the x,y,-plane when the vane lies entirely inside 
the Mach cone is shown in Fig. 2. Transferring to the f-plane, the tips of the vane 
x,=0, y,= +c are represented by t= +ia respectively where 


The boundary conditions could be laid down directly in the ¢-plane; they are: — 


a: (20) 


real axis, for 1<0 and say 
imaginary axis, for |t|<a. 
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Fig. 5. Fig. 6. 


These boundary conditions show clearly the convenience of transforming the 
cut t-plane to a cut z-plane by the transformation 


To make the transformation from the f-plane to the z-plane one-to-one, z is 
taken to be real and positive when ¢ is real and positive, and by following the 
argument of z 


arg z=4 [arg (t+ ia) + arg (t—ia)], 


gives for the upper half-planes : — 


(i) The real axis in the t-plane is transformed to those parts of the real axis 
in the z-plane for which | z|>a. 


(ii) For the vane OA (Fig. 6), the side # (t)= +0 is mapped on to that part 
of the real axis for which 0 <z<a and the side #(t)= —0 is mapped 


on to —-a<z<0 on the real axis. The point A is mapped on the origin 
in the z-plane. 


(iii) The imaginary axis with | t| >a is mapped on the whole imaginary axis 
in the z-plane. 


The lower half is transformed similarly. The correspondence between the 
points on the real axis is shown in Figs. 5 and 6. 


The boundary conditions in the z-plane are all given on the real axis and are 
as follows :— 


v= —-Wa for say 
and for z<—a, 
v=0 for z, >z>(1+a*}i=b say, 
u=0 for —a<iz<a, 
w-0 for a<iz<b. 


Since u,v, w are the real parts of functions of G,(t), G.(1), G, (0) respectively, 
they are also the real parts of functions of z, H,(z). H.(z), H(z) respectively, so that 
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u+iu, =H, (2), 
w+iw, =H, (z). 


Hence, from equations (14) and (22), 


H,’ (2): Hy’ (z)=1: : —tan p (25) 


The boundary conditions (23) lead to the following conditions on H,’ (z):— 


(z)} —0 for bh<z and z<~-a, | 
(26) 
J (H,'(2))-0 for ~-a<z<b. | 
The conditions (23) and (26) show that H..’(z) must have the form 
m+ — 


(z—2Z,) 


where m, n, | are integers and H(z) is an integral function which is real when 
z is real. 


The following restrictions must be imposed in order that the solution may 
correspond to a physically possible flow. 


(i) In the physical plane v= — Wz at x,=0, y, = +1 and in a neighbourhood 
of each of these two points, and therefore as in Ref. 1, H,’(z) must be 
O(z-*) at most as | z|—> ©. 


(ii) u,v, w must be finite at z=b and z=-—a. As H,’(z) contains the factor 
(z—b)"*+, H,’(z) and H,’ (z) contain the factor (z— 5)" and so n= 0; also 
it can be proved that m>-1. 


(iii) At the subsonic leading edges of the vane, the velocity normal to the plane 
of the vane must vary as r,! as r,—>0, where r, is the normal distance 
from the edge in the plane of the vane. Therefore in the f-plane, 
G, (t) ~ K (a+ 1’)! as t—> ia or — ia, the two points corresponding to the 
two edges of the vane (K is 4 constant). In the z-plane H, (z) ~ K/z and 
therefore H,’(z)~-—K/z? as z—>0; hence /=0. 


(iv) If condition (i) is to be satisfied, then H(z) must be a polynomial of 
degree p say. The boundary conditions show that there must be no 
residue at the double pole z=0; hence p> 1. 


Thus from these four conditions 


m+n+p<0, 
m= -1l, 
n>0, 
p=x-l, 
262 


4) 


6) 


Ly 


VANES ON SUPERSONIC WINGS 


and the unique values of m, n, p satisfying these inequalities are 


m=-1, n=0, p= 1. 
(b 
Therefore H,’(z)=A (z+B), . (28) 


where A and B are real numbers. The value of B is obtained by equating the residue 
at z=0 to zero and that of A by equating the residue at z=z, to iW2/x, whence 


2abz, 
= (a+b)—2ab’ (29) 
_ Wafa+z,\! (a+b)—2ab 
( —b a+b (30) 
From equation (25) 
z+B z+a\! 
, B 


Hence in the upper half-plane 


(Gra) 2 (atz 


H, (z)= iWa — py a+b\z,—b z 


where the square roots are positive when their arguments are real and positive, and 
the inverse tangents are between — 7/2 and ~/2 for real arguments. 


The velocity components are obtained by taking the real parts of these 
expressions for the H’s after substituting for z in terms of ft, and thence in terms of 
X,, ¥,- For general points inside the Mach cone this is a very tedious process, but 
for points on the wing and the vane, that is for points on the real axis in the z-plane, 
the results follow immediately from equations (33) to (35). 
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Thus on the upper wing surface w is given by the right hand side of equation 
(35), which is purely real in the interval —b < z<—a, where now 


c 1 


(a’ + =(a* + sec? o,)}, a b 


Case 2. 1<c<cosec u, 


In this case the edges of the vane are outside the Mach cone. Therefore the 
fields in the two upper quadrants are independent, and in the f-plane they will be 
separated by the whole imaginary axis. The two fields have then to be solved 
separately. 


For the left-hand quadrant in the ¢-plane the boundary conditions are 
v= —Wa for —seco,<r<0, 


and from (14) (36) 
v= —Wa(l+coto,tano,) for t< —secv,. 


The condition u=0 on the vane in the x,y,-plane becomes # {G,’ (t)} =0 on the 
imaginary axis. Hence on the imaginary axis .7{G,’ (t)} =0 and from equation (19) 
R{G,/ (t)}=0. Thus the boundary conditions for G,’ (t) are 


(t)}=0 on the real axis for 1<0 say), 
) 


R{G,'(t)}}=0 on the imaginary axis. 
The condition that G,’(t) is O(t-*) at most as t—> 00 still applies, and the dis- 
continuity at t= —?, in the value of v along the real axis shows that G,’(#) has a 


pole there with residue (iWa/x)coto,tanc,; G,’(t) has no other singularity in 
the upper left-hand quadrant of the t-plane or on its boundary. 


The solution for G,’(t) and G,(¢) with these conditions is unique and is 


= 1 


From equation (19) 


G,’ (t)= cote, tan (39) 
Hence in the upper left-hand quadrant: 
_ Wa = tht} 
G, (t= cotc, cos cos (41) 
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iWz t—t 
G, -- —— cot, tanv, log - W2z(14-cote, tane,), (42) 
t+t, 
Wx { _, _, 1+t,t) 
G, ()= tan» cot«, seca, COS + cos * 


+ Watan »coseco, (43) 
The real part of these functions could be obtained easily. 


The value of w on y,= +0, — L<x,<0, i.e. on the upper surface of the 
wing is 
w= — tan sec«r, COs — + cos 
+ 
+ Wa tan » cosec (1 +cos«, sec«r,). (44) 
In the upper right-hand quadrant, the boundary conditions on G,’ (t) are 


JS{G,'(t)}=0 for O<t<l, | 


for 1<t, 
when ¢ is real, and (45) 


=0 for |t|>0, 
when ¢ is imaginary. 


With these boundary conditions, following the method used previously, it is easy 
to find that 


Integrating, in the upper right-hand quadrant 


t—t, 
t+t, 


} —Wztan u(t, —t,) coto,. 


iWa, f 
G, (= tanpcote, log —t, log 
(51) 
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Case 3. c>cosec a, 
In a similar manner to Case 2, for the upper left-hand quadrant 


2t, 

G, (= P—t,2’ (53) 
Ww 2t, t 

G/()= —tanp (54) 


Integration gives, over the surface of the body, 


Wa 


2Wz tan cosecc,. 
to +X, } + 


(55) 


For the upper right-hand quadrant u vanishes on the boundary; hence G, (1), 
G. (t), G; (t) vanish everywhere in this quadrant. 


w= tan p { cos-* 


Having obtained the components of the perturbation velocity in the upper half- 
plane, outside and inside the Mach cone, the corresponding values of these com- 
ponents for the lower half-plane are obtained for the three cases from the fact that 
v is even, while u, w are odd, in y,. 


5. Lift Coefficient 


The lift force in the direction of the positive y-axis is obtained by integrating 
the perturbation pressure over the upper and lower sides of the wing. When the 
wing and the vane have straight trailing edges normal to the direction of the 
undisturbed stream, the lift coefficient (referred to area and 4pW7’) for the part of 
the wing inside the Mach cone of the vertex is 


a, . ‘ (56) 


2{ 


-1 


where P is the perturbation pressure on the upper surface given in equation (2). 


Therefore | wdx,. 


-1 
On the wing, x,=t, w=G,(t), (-1<t<0), 


0 


and hence C= wl G, (t) dt. 


-1 


266 


(52) 


(53) 


(54) 


(55) 


1alf- 
om- 
that 


ting 
the 
the 
t of 


(56) 


(2). 


VANES ON SUPERSONIC WINGS 


Integrating by parts, 


0 


C,=42 cosec o, tan — wl G,’ (t)tdt. . (57) 


-1 


For the first case, c<1, transform to the z-plane. From equation (22), 
z= —(a*+12’)! on the wing, so that -b<z< —a. Hence 


C,= 42 cosec o, tan u+ H,,’ dz. . (58) 


An expression for H,’ (z) is obtained for this case in equation (32). Substituting and 
then performing the integration, 


cosec o, tan { 1- [ 2. (0) ~ 


(59) 


When the leading edge of the wing is perpendicular to the stream direction, z,=00 
c=" /2, 


2 
C,=42 tan [1- ] 


which is the value obtained by Stewartson®’ who solved this problem by using a 
Laplace transform method. 


For Case 2, 1 <c<cosecr,, substituting for G,’ (t) from equation (40) in (57), 


= 4zcosec tan [ cos o, tane (ata 
and c=cosec 
c 
Therefore C,=42coseco, tan [ 1+cos¢,, . 


For Case 3, c>coseco,, C;, can be obtained as in Case 2, giving the result 
C.=42acosec o, tan » [1 + tan +, (cosec — 1)]. (62) 


If C,’=4a cosec r, tan » is the two-dimensional value of the lift coefficient, the 
variation of (C;—C,’)/C,’ as a function of c is calculated for three different wings, 
having values of sweepback of 0, tan~' (tan u/2) and uy, i.e. for sec! 2, 
sec-' 1). The results are shown in Figs. 7 and 8. 
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Fig. 7. 


It will be seen from these figures that the effect of the vane is to reduce the 
deficit of lift otherwise caused by the presence of the wing tip. With zero sweepback 
the value of the lift coefficient tends to the two-dimensional (Ackeret) value as the 
vertex semi-angle of the vane increases towards the Mach angle, i.e. as c increases 
from 0 to 1. When the wing is swept, the effect of the vane is similar. Its effect is 
to restore the lift to the two-dimensional value for a vane semi-angle, y=y, < /, 
depending on the angle of sweep 8; so that for 0<y</,, the lift is less than the 
two-dimensional value and for y > y, the lift is greater than this value. It reaches 
its maximum when the vane edge touches the shock wave produced by the leading 
edges of the wing. 


This result was encouraging, and it was decided to go further and determine 
the drag and to discuss how the overall performance was affected by the presence 
of the vane. 


6. The Symmetrical Part 


This part corresponds to the wing and the vane of small finite thickness, but at 
zero incidence; the drag exists because of the thickness. The form of wing and vane 
for this case has already been described and a sketch of it is shown in Fig. 1. The 
thicknesses of the wing and the vane are represented by the angles 4 and < which 
their surfaces make with the z-axis. The solution is given here only for the field 
inside the Mach cone when the edges of the vane lie entirely inside the Mach cone. 
The other cases can be treated in the same way as the corresponding anti-symmetrical 
cases if they are needed. 


To the degree of approximation assumed, the boundary conditions may be 
applied on the planes of symmetry of the wing and the vane instead of on the 
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Fig. 8. 


actual surfaces. Then, by neglecting quantities of the second and higher orders, the 


boundary condition on the upper surface of the wing y,= +0 is 


v= Ws, 
and on the lower surface of the wing y,= —0 is 
v= —Ws, 


while on the right-hand side of the vane x,= +0, the boundary condition is 


u=We, 
and that on the left-hand side of the vane x,= —0 is 
u=—We, 


where W is the velocity of the main stream in the z-direction. 


The velocity components on the Mach cone are similar to those obtained in the 
lifting case in equation (11) but with z= —3 in region 1 and a=8 in region 2. 


On transforming into the t-plane defined in equation (17) and then to the z-plane 
defined in equation (22), the components u, v, w of the perturbation velocity are the 
real parts of the functions H, (z), H, (z), H, (z) of the complex variable z. The bound- 


ary conditions could easily be transformed to the z-plane; they are 


u= F&F {H, (z)}=We for 0<z<a, 
=—-W: for -a<z<0, 

v= #{H,(z)}=Ws for z,<z and for z<—a, 
=Q for a<iz<%, 


269 


0-6 
| SEC 
a 
0-4 
Z 
/ 
/ 
c 
the 
ack 
the 
ses 
t is 
the 
hes 
ing 
ine 
1ce 
at 
ine 
‘he 
ne. 
. (63) 
be 
he 
|__| 


M. A. 


GORGUI 


and hence from equation (28), the boundary conditions on H,’ (z) are 


J {H,/(2}=0 for -a<z<a, 


(64) 
R {H,'(z)}=0 for z<—a and z >a. | 

From these boundary conditions H,’(z) must have a double pole of residue 2iW: /= 

at the origin, since it corresponds to the subsonic leading edge of the vane, another 

pole with residue iW4/« at z=z,, and a singularity at z= +a of the form (a? —z*)"*! 

where n=>-—1. Therefore H,’(z) must have the form 


(a? 


H(z) . (65) 


where H (z) is an integral function which is real when z is real. Since H,’(z) must 
be of O(z-*) as z—> ©, H(z) must be a polynomial of degree p say. Therefore 
we have 

2n+p<0, 


n>-l, 


p21. 


Then n= —1, and either p=1 or p=2. In either case H,’ (z) can be put in the form 


(66) 
and from equation (27) 
H,’(z)=-i (67) 
H,’ (z)= —tanp (S ~ (68) 
In these expressions 
2We ab, (z,? —a?}. ° # (69) 


The value of C is found after integrating H,(z), from the condition that 
u=& {H, (z)} vanishes for b<z<z. 


The expressions (66), (67) and (68) are now integrated. H,(z) and H, (z) can 
be expressed as elementary functions. H,’(z) can be written in the form 
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On integration, the first and second terms give elliptic integrals, and the remain- 
ing terms give elementary’ functions. 


For the first and second terms, let 


1 a 
k= k’=(1 —k?) : . (71 


x is defined by means of Jacobian elliptic functions 
dn (x, k)=z/b, sn (x, k)=(b? — cn (x, . (72) 


The elliptic integral of the second kind 


E (= | dn? (6,k) a9, 


and the elliptic integral of the third kind 


_sn? 
b? 


If x=K, the complete elliptic integral of the first kind to modulus k, E (x) and F (x) 
will tend to the corresponding complete integrals E and F. The values z= —0©o, 0, 
a, b correspond with x=iK’, K+ iK’, K, 0, where K’ is the complete elliptic integral 
of the first kind to modulus k’. With all these substitutions, the integral for H,’ (z) 
on the upper half-plane that vanishes at z=a is 


C 2 g?)t — 72) WS 

H.@= W-B] - FW-P- 
Ws 2? ( We b(a?—2)t+a(b?— 2)! 
3) (4-5) log +2We. 
(75) 


On the real axis for b<z<z,, x, E(x), Fx) are all purely imaginary, and in 
order that .%{H, (z)} =2Wée in this range, 


C_ Wi *F 

Therefore H,(z)= Z bE | 


2iW log 


+2W:. 
(77) 
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The integrals for H(z) and H, (z) on the upper half-plane are 


(78) 
C 2Wd/( z,?-a@ _, 
2Wea b+(b?— 


The velocity components may now be obtained by taking the real parts of these 
expressions after substituting for z in terms of t and thence in terms of x,, y,. The 
corresponding values for the lower half-plane are obtained from the fact that in this 
symmetrical case v is odd, while u and w are even in y. 


7. The Drag Force 


The drag force in the z-direction is due to the perturbation pressure acting 
normally over the surfaces of the wing and the vane, and to the suction force normal 
to the leading edges of the vane. The first two parts, the drag over the wing and 
the drag over the vane, are obtained by integrating the z-component of the force due 
to the perturbation pressure over their respective surfaces. The third part is obtained 
by integrating the z-component of the suction force over the two leading edges 
of the vane. 


The drag over the wing 

Let P be the perturbation pressure defined in equation (2). The drag over that 
part of the wing inside the Mach cone which lies between the origin and the plane 
z=cot pw, i.e. the drag over an area }cot» of the wing, is 


{ (at }. 
-1 


Integrating by parts 


0 


D,=-—pcotuW3F { Ke 


Changing the variable of integration from f’to z, from equation (24), z= —(a’+ t’) 
on the wing, so that the range of integration is -b<z<-—a. Hence 


D, =p cot { [ (z)(z?- a) | H,’ (2) (z2— dz. \ (80) 
The expressions for H,’(z) or H,(z) are all real in —b<z<~—a. 
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The drag on the wing is composed of two parts; the first part is due to the 
incidence of the wing, and for this H,’ (z) and H, (z) are given by equations (32) and 
(35) respectively. The evaluation of (80) in these circumstances was carried out 
when the lift on the wing was calculated in equation (59), and hence 


D,’ =pW?2? cosec o, {1 (Be ( (81) 


The second part is due to the thickness of the wing and the vane; and here 
H,'(z) and H,(z) are given in equations (68) and (79) respectively. Therefore 


f WS (z,?—a*) 
u| W8cosec tan tan ((S + )( >) dz 


2 
b? 


= pW28? [ coseco,+ gr (K-E)+ (K—F)+ 


The total drag on the wing is D,’+D,”. With no sweepback z,—> 00, 
«,—>7/2. From the definition for F and K given before it is easy to prove that 


Lim F=K, 
Therefore cu K 
E 
and 
D,=pWs [ 1+ ( E E)|+ow a Aq) 
(83) 


The drag over the vane 
In this case the pressure is to be integrated over the two sides of the vane, from 


the vertex at the origin to the plane z=cot u, i.e. over an area ccotu. The drag 
will be 


c 


D,=— pe cotuW | ow 


+w 


a 


) dy. 


+0 
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On the vane, if 7’ is the imaginary part of t, from equation (17) 


dt’ 
Therefore —pWe cot | (Wre (ty=40 + Wre (t)=-0) 


On transforming to the z-plane, z=(a’ —1’*)!, 0<( z<a on the right-hand side 
of the vane, and z= — (a? -—?’*)', —a<.z<0 on the left-hand side. The drag will be 


6 
D.= — cot | + 


(b? (a? — : (84) 


From the expression (79) for H,(z) it is clear that it is real in the interval 
—a<z<a. Therefore 


pWeC 
+ f (b? — (a? — 

0 0 

2pW?e7a | a+ — zdz = 
— log (b? — (a? — 


Integrate the first term by parts, then divide the integral part into its even and odd 
parts. The even integrals in D. will vanish on integration. Calculating the odd 
integrals gives 


W726 


Zo b 
a 
= | —tan~ 'a} + 
4pW?22q 
a J 


The last term is the drag over the vane in the absence of the wing, a result 
obtained by Puckett“). The remaining two terms represents the effect of the thick- 
ness of the wing. 
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The drag due to the suction force on the leading edges 
of the vane 
At the subsonic leading edges of the vane, the velocity becomes infinite in the 
linear theory both in the anti-symmetrical and symmetrical part of the field. If at 
any point (z, y) on the upper leading edge of the vane, r, is the normal distance from 


Fig. 9. 


this edge in the plane of the vane (see Fig. 9) the velocity normal to that plane 
~zr,~? as r,—>0, ie. u~:,r,~? (where <, is a constant proved later to depend 
on y); there is a suction force on the leading edge and in the plane of the vane, whose 
magnitude per unit length is“? 


The drag force due to this suction force over a length cosec y of the upper edge 
of the vane is 
Tpe,” 
sin? 
In conical flows the components of the perturbation velocity u, v, w are constant 
along straight lines emerging from the vertex of the Mach cone. Let ¢,’ be the 
value of <, at the upper edge of the vane in the plane z=coty, i.e. at the point 
(cot n, c) in Fig. 10, so that u~<,’/r,t as r,,—>0. Therefore along any straight 
line, A say, near the upper edge of the vane, 


 & 
Hence (24) (2) ‘ . (87) 


y 
gi 
(cot 
oS 
Se 
(1 — M? sin? 
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Therefore the drag force due to the suction force along the upper edge of the 
vane is 


~ 2(1—M? sin? 


A similar expression is obtained for the lower edge, and the total drag due to the 
suction force becomes 


TPC (¢,7+ 
2(1 — M? sin? y)’ 


D,=—- (88) 


where «,’ is the corresponding value of <,’ at the lower edge of the vane. 


To determine <,’, 
H,’ (z)~ as z—>0 along the imaginary axis. 


Then in the t-plane, 


G’O~ 


as t—>ia along the imaginary axis. 


In transforming to the x,y,-plane 


dy, (y,?-—c*)/? (l—y,?}! 


as y,—>c, 


(1 1 
_,a-e) 1 9 
Hence u (89) 
Comparing equations (86) and (89), 
e,’= — cost y (1 —c?)(2c)-*A. (90) 


But from equations (32) and (68) the value of A corresponding to the upper edge of 
the vane is 


Wa C 
a+b \z,-b ab’ 


and that for the lower edge is 


Wa | Cc 
a+b \z,-b ab* 
2 
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Therefore the drag force becomes 


Since c=tan y cot u, M=cosec hence 
ay] Cc 2 (ab)! Z +a ey 


The total drag is now obtained by adding together the expressions obtained for 
D,, D., D,. For no sweepback this gives 


(SE -E)| +1)+ 


4c sin” ] | - 
+ — sin-'¢ 98° + pW*a?) 1 2 | 


[Ke 4a? 


14 


When 2 is very small so that z* can be neglected, 


c(2+c)« , cloge | 
| sin le: (93) 


The variation of the drag given in this last equation as a function of 
c=tany/tan can be seen in Fig. 10. 
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8. Performance 


It is now possible to discuss how the performance is affected by the presence of 
the vane. Let L be the lift force at incidence « and D the drag force at zero 
incidence; the trailing edges of both the wing and the vane are straight lines lying 
in the plane z=cot . Then the quantity 


=> , ‘ ; (94) 
( 
is a measure of the performance of the wing-vane system. 

Consider the only available case, that is when the vane lies entirely inside the 
Mach cone. The other cases when the vane edges are outside the Mach cone are 


not important, since in these cases the drag increases without any considerable 
increase in lift. 


7 is obtained by substituting for L and D from the values already obtained. To 
avoid unnecessary complications the case of no sweep (z,=0, +,=7/2) is 
considered. 1 is then given by 


(95) 


The variation of 18*/a as a function of c, given by equation (95), is shown in 
Fig. 11. The calculations were made for wings carrying vanes of different thickness 
in a flow for which the Mach number was /72. 
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It appears, then, that with no sweep, an increase in the performance is made 
possible by using wing tip vanes, provided that : < 4; for each value of < there is an 
optimum semi-angle for the vane, this angle increasing as « decreases. The maxi- 
mum performance is obtained when <=0, which clearly must be the case, the index 
of merit » then being about 1-1 2/8’ and the semi-angle of the vane being about 25° 
when the Mach angle is 45°; with no vane »=0-735 2/8", so that the performance 
of that part of the wing-vane system which lies inside the Mach cone is increased by 
a factor of 1:5. In the more practical case when <=0-238, the optimum value of » is 
0-925 a/3? for a semi-angle of about 18° when the Mach angle is 45°, the perform- 
ance in this case being better than that with no vane by a factor of 1:26. Or when 
¢=0°55, the optimum semi-angle is about 65° when the Mach angle is 45°, and 
the factor of improvement is 1-07. 


It is clear that an increase in performance also occurs when the wing has a 
swept leading edge, and the indication of Fig. 7 is that this may be greater than in 
the unswept case. 
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Estimation of the Effects of Distortion on 
the Longitudinal Stability of Swept Wing 
Aircraft at High Speeds 
(Sub-Critical Mach Numbers) 


B. S. CAMPION, B.Sc.(Eng.), D.C.Ae. 
(College of Aeronautics, Cranfield) 


Summary: The effects of distortion on the longitudinal stability of swept wing 
aircraft at high speeds (sub-critical Mach numbers) are considered on a quasi-static 
basis. The method is based on the theory of Gates and Lyon") but involves some 
extension of this theory. 


The treatment of wing distortion is considered in some detail and the effects of 
built-in twist and camber and wing weight are included, using the so-called superposition 
method.) The application of the analysis of Lyon and Ripley®) for investigating 
fuselage, tail and control circuit distortion is suggested, but means of modifying and 
simplifying this procedure where desirable are put forward. 


The analysis is illustrated by means of a simple example of a swept wing fighter 
aircraft for which wing, fuselage and tail distortion effects are considered, and the 
results are discussed with reference to the individual and combined distortion effects, 
as well as the effect of compressibility. 


Introduction 


The longitudinal stability of an aircraft is usually considered in terms of “ static 
stability” (measured by the static margins), “ manoeuvrability” (measured by the 
manoeuvre margins) and general dynamic stability. If the stability derivatives are 
very little modified by frequency effects (i.e. they approximate closely to the quasi- 
Static values), then the static and manoeuvre margins as normally defined can be 
related to the coefficients E, and C, in the “stability quartic,” given by the usual 
small displacement theory. The values of the static and manoeuvre margins then 
largely determine the characteristics of the phugoid and short period oscillations. 


When structural distortion effects are introduced, it is again necessary to 
consider whether the motion of the aircraft (in the long and short period oscillations) 
occurs under quasi-static conditions, or whether the dynamics of the separate aircraft 
components should be considered, introducing additional degrees of freedom 
corresponding to oscillations of individual components. This question is discussed 
in Ref. 4, where it is pointed out that if the frequency of the short period oscillation 
is of the same order as the lowest natural frequency of any component (e.g. the 
wing) the simple quasi-static approach is suspect. Once the quasi-static approach 
is abandoned, however, the treatment of dynamic stability when distortion effects are 
included becomes very difficult. The general treatment of the dynamic stability of 
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a flexible aircraft can be similar to that used in flutter problems, although the 
difficulties are enhanced by the fact that coupled oscillations of wing, fuselage, tail, 
and so on, are combined with the overall “rigid body” motion of the aircraft. An 
attempt to formulate the equations governing the motion of an aircraft with flexible 
fuselage and wings is made in Ref. 5, and the problem is considered briefly in Ref. 4. 


In this report the “ quasi-static” approach only is considered, as in Refs. 1 and 
2. By this method the equations of motion for a rigid aircraft are used but the 
values of the aerodynamic derivatives are modified to include distortion effects. 


The basic theory is, in essentials, that of Refs. 1 and 2, with certain modifica- 
tions and extensions which it is believed will permit of a more logical treatment of 
the effects of wing distortion in cases where there is built-in twist or camber, and 
the treatment also permits the ready inclusion of the effects of aircraft weight. The 
conditions for the balance of aerodynamic, elastic and inertia forces are obtained by 
the superposition method of Ref. 3, which it is considered has many advantages over 
other methods. 


The analysis is illustrated by means of an example of a high-speed swept wing 
fighter, and the results of the analysis are discussed in detail. This example presents 
a number of features of general interest. 


For compressibility effects on a rigid aircraft the important parameter is the 
Mach number, so that the variation of true air speed must be considered. When 
the aircraft is flexible, however, the distortions produced by aerodynamic loading 
are dependent on the equivalent air speed. In general, therefore, the two parameters 
M and q=4pV* must be considered. In this report the suffix M indicates that a 
derivative is taken at constant Mach number (e.g. (@C;,/0)y), and in such cases it is 
also implied that g is constant. 


The corrections for variations of inertia loading due to normal acceleration, 
introduced in Section 4, are “quasi-static,” i.e. it is assumed that the normal 
accelerations of all parts of the aircraft are the same as that of the c.g: and that the 
structure is always in equilibrium under the applied aerodynamic and inertia loading. 
This assumption is similar to the assumption that frequency effects on the aero- 
dynamic derivatives, and so on, may be neglected and will similarly become invalid 
when the short period frequency is high. 


Notation 
A _ wing lift curve slope with compressibility and distortion effects 
included 
Ay rigid wing lift curve slope when compressibility effects are 
included 


A° value of A for condition n=0 (i.e. neglecting inertia loading) 
An (OC inw/ 

As, Bs,Cs; superposition coefficients 
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AvhwA, | tail lift coefficient derivatives’ when compressibility and all 
ams distortion effects are included 
Aj1, Byr = lift and hinge moment derivatives corresponding to true tail- 
setting angle 
0A, 
’ 0M defined in Section 4.3 
etc.- 


B,,B,,B,,B, tail hinge moment derivatives when compressibility and all 
distortion effects are included 


C standard mean chord 


Ci, overall lift coefficient based on wing area in steady level flight, 
= W/(4pV?S) 


Ciw wing lift coefficient 


Crews CitwsCunw Wing lift coefficient increments due to built-in camber and twist 


and wing weight 
Crow wing lift coefficient due to incidence only 
Ciwe Crew Crtw 


Ci» Wing pitching moment coefficient about rigid wing mean aero- 
dynamic centre 


Crews Cmiws Cw Wing pitching moment coefficients corresponding to Crew, Criw. 


Law 


(Cmo)wing | Wing pitching moment coefficient about flexible wing mean aero- 
dynamic centre 


C+ fuselage pitching moment coefficient (based on wing area) 
Cy = W/(GpV’S), “resultant vertical force coefficient ” 


F= 4%) 


Fr. 7 ) 
F = 


H,,,H»’ — stick fixed, stick free manoeuvre margins 


H,,H,° value of A corresponding to position of (flexible) wing mean 
aerodynamic centre of lift due to incidence only—affix 0 refers 
to n=0. 


H," “inertial aerodynamic centre” position 
H,, value of H, for rigid wing 
K,,K,’ _ stick fixed, stick free static margins 


| distance from tail aerodynamic centre to wing aerodynamic 
centre (corresponding to H,°) 


ly distance from tail aerodynamic centre to c.g. of aircraft 
n defined by (n—1)g=normal acceleration 
q = pV’; also angular velocity in pitch 
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= 
Sc tail 
+ volume 
Paw Vac. | coefficients 
1+F 
: 1+F 


a wing incidence measured at root 
Gq, %qp, tailplane root incidences corresponding to ny, nr, 
8 tab angle corresponding to movement of pilot’s trimmer control 
€ mean downwash angle at tailplane 
€, mean downwash angle at tailplane due to wing root incidence 
€, mean downwash angle at tailplane corresponding to Ciw, 
é, mean downwash angle at tailplane due to distortion of wing 
under inertia loading 
= €./(A°2) 
= €,/(A,n) 
» elevator angle corresponding to movement of control column 


nro _ tail setting angle corresponding to setting on ground (if “ fixed ”) 
or to movement of pilot’s control (if “ variable”) 


nr _ true (root) tail-setting angle 


2. Wing Distortion 


When considering longitudinal stability it is usual to assume that the ailerons 
will remain in the neutral position when the wing distorts. Distortion of the main 
wing structure only is then considered. For most swept wings the loss of incidence 
due to upward bending is greater than the increase due to twist (referred to an ideal 
straight flexural axis). Thus wing distortion produces a net loss of lift curve slope, 
accompanied by forward movement of the aerodynamic centre (see Section 2.2). 


When the quasi-static approach is used in obtaining stability criteria and 
derivatives it is assumed that the structure is always in equilibrium under the applied 
loading. The problem of treating wing distortion is then basically that of solving 
for the “ aeroelastic equilibrium ” of a flexible lifting surface. 


2.1. Solution of the aeroelastic equilibrium problem for a flexible 
lifting surface 
The problem of calculating the aerodynamic characteristics of a flexible lifting 
surface is one of some difficulty, due to the fact that distortion under load produces 
a change of load. Mathematically the problem takes the form of the solution of an 
integral equation. 
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It is possible to solve the problem (using strip theory) by successive approxima- 
tion (see Ref. 4) or by the method of semi-rigid representation “:®. These methods 
are well known and widely used. They suffer from the disadvantage that, since strip 
theory must be used, induced aerodynamic effects due to the distortion itself are 
neglected; also with the semi-rigid method, as normally applied using one mode, 
the accuracy is limited by the closeness with which the mode represents the distortion. 
However, the semi-rigid method is not in principle limited to one mode and the 
“superposition” method which is put forward in Ref. 3*, and which has been adopted 
here as being in some respects superior to previous techniques, can strictly be 
regarded as a variant of the semi-rigid representation technique. The following is 
a brief description of it. The usual assumption of linearity between loading and 
incidence is fundamental to the method. 


At any point on the span of an elastic wing in equilibrium under aerodynamic 
load, the final geometric incidence zy will be given by 
2+ 
where geometric incidence of undistorted wing (“ initial” incidence) 


%, change of incidence due to distortion. 


If zp is known, the final loading is known and 2, may be calculated, giving 2. 
Thus the problem is easily solved “ backwards.” If a number of arbitrary “final” 
incidence distributions zp,, z7., and so on, are chosen and the corresponding 
“initial” incidence distributions 2;,, z;,, and so on, are so obtained, then any 
given initial incidence distributions may be approximately represented by a linear 
combination of the arbitrary initial distributions. Thus 


Asay, + + Caz, + 


In practice three or four such terms may be sufficient. 


The coefficients As, Bs, Cs, and so on, are functions of the parameter gAx, where 
Ax is the rigid wing lift curve slope allowing for compressibility effects. 


But Ag2z,= As@p, — 
and therefore Ap — Agap, + + + 
— (Asp, + +...) 
so that Agap, + 
and hence for the given initial incidence distribution, «;, the final incidence distribu- 
tion 2p, and hence the final aerodynamic characteristics, can be obtained. By means 


of this method the aerodynamic and structural problems are separated and may be 
considered independently. 


*The superposition method has also been developed independently at Handley Page Ltd., but 
this work remains unpublished, 
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This method is likely to yield more accurate results for many problems than the 
other methods mentioned previously, due mainly to the fact that induced aero- 
dynamic effects due to distortion are readily included and no sweeping assumptions 
need be made about the mode of distortion. Once the calculations have been 
completed for the arbitrary cases chosen, equilibrium conditions for any combination 
of gq and M are quite easily obtained. As with other methods, some difficulty is 
encountered if the form of the aerodynamic loading varies appreciably with Mach 
number (see Section 5.2). 


2.2. Effects of wing distortion on the wing lift and pitching 
moment contributions 


2.2.1. WING WITH ZERO BUILT-IN TWIST AND CAMBER; DISTORTIONS DUE 
TO WING WEIGHT NEGLECTED 


If the incidence of the wing < is defined as the angle between the chord line of 
the wing and the direction of flight, measured at the wing root, then 
a M 
where Ciaw= Wing lift coefficient of flexible wing due to (root) incidence x 
. 
A= wing lift curve slope (following Ref. 1)= et) ‘ 
M 


It is possible to find (@C,aw/02)y using the superposition method as already described 
(see also Ref. 3). 


Similarly, taking moments about the rigid wing aerodynamic centre, 


where Craw is the pitching moment coefficient corresponding to Ciaw and may also 
be found by the superposition method. For the “rigid wing” Cnw as defined is of 
course zero, but with the flexible wing there is a pitching moment about the rigid 
wing mean aerodynamic centre, H,x, which is proportional to incidence, i.e. there 
is a movement of the wing aerodynamic centre given by 


4H, = 0a 0a 


(3) 


giving the mean aerodynamic centre of the flexible wing 


H,=H,2+ 
285 


a- 
rip 
are 
ade, 
on. 
he 
ed 
be 
is 
nd 
ic 
” 
y 
r 
t 


B. S. CAMPION 


The wing pitching moment coefficient about the aircraft c.g. is then 


Caw a= Aa (h H,) (4) 


and the wing pitching moment coefficient about the new mean aerodynamic centre 
H,, is zero, 


i.e. = 9. 


2.2.2. WING WITH BUILT-IN TWIST AND CAMBER; DISTORTIONS DUE TO WING 
WEIGHT NEGLECTED 


When the undistorted wing has twist or camber it is possible to consider the 
effects of incidence, twist and camber separately. Applying the principle of super- 
position as in Ref. 3, 


OCraw 
( ) a+ + Citw 


(5) 


a=root incidence as defined in Section 2.2.1 


Ciow= lift coefficient due to root incidence z on wing with zero built-in 
twist and no built-in camber 


Cicw= lift coefficient on wing with zero root incidence and zero built-in 
twist, but with built-in camber 


Citw= lift coefficient on wing with zero root incidence and zero built-in 
camber but with built-in twist 


Case Ci we + Cy tw 


(as before). 


Crew and Cyt may be found using the superposition method and are functions of 
speed and Mach number. 


Thus cm) = = Cun) =A as in Section 2.2.1. 
M 


Also, the pitching moment coefficient about the rigid wing mean aerodynamic 
centre is 


0a 


where Crows Crows Cmtw are the pitching moment coefficients corresponding to Crow. 
Crews Crtw, and hence 
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ie. AH,= / as in Section 2.2.1. 
Then Cawos.= %+ + (h— Ayr) Ciws 
where Cuwo= Carew + 
i.e. A2Q4H, + +(h— Hon) (A2+Crwo) 
= (Cino)wing + A@ (h—H,), ‘ (6) 


where (Cro) wine Cue + (h Hi») 
and H,=H,p+ Ay. 


This expression for Crwe,¢. is different from that normally used in that Cry Az. 
It will be seen, however, that it is convenient to express quantities in the form 
X+Ya, where X and Y are functions of speed and Mach number but not 
of incidence. 


2.2.3. WING WITH BUILT-IN TWIST AND CAMBER WHEN DISTORTIONS DUE TO WING 
WEIGHT ARE INCLUDED 


Following the method adopted in Section 2.2.2, 


OC nw 
where n—1= (1/g)x normal acceleration of aircraft (positive upwards) 


i.e. n= C,ApV?S/W, in flight for which the inclination of the flight path 
to the horizontal is small 


Cinw= lift coefficient on wing with zero root incidence, zero built-in twist 
and no built-in camber, due to deflection of the wing under its 
effective weight. 


The affix 0 denotes the condition n=0, so that the effects of wing weight are entirely 
contained in the last term. 


As with Cyw,, and so on, 0C;jn./0n may be found using the superposition 
method. 


Similarly, taking moments about the rigid wing aerodynamic centre, 


awe 
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and about the c.g., 


n 
ec, nw 
+ (h Hi») (A°%z (h H,, ) 
n M 
0 OC aw y 
An= ( on on on 


Ayr 


In practice the affix 0 could be dropped from A°, Cywo°s Cmwo® aNd (Cno)wing? Without 
confusion occurring. 


The effect on existing longitudinal stability theory of using these more compli- 
cated expressions for Cin, Cmweg. is considered in Section 4. 


2.3. Effect of wing distortion on the downwash at the tail 


The incidence of any point on the rigid tailplane is given by 
where ¢,=local downwash angle. 


Hence, following the method used in dealing with wing distortion, 
ACiz=A ya nr)—A 1e=R 


due to tailplane incidence, 


where Aja= 02, 
%,= root incidence of tailplane with zero wing downwash 
Nr 
A je= OCi7/ 
&,=downwash angle at tailplane root. 


Values of A,. and A,. might then be found using the superposition method*. 
provided that the spanwise downwash distribution were known. 


*This is strictly true only for an all-moving tailplane (see Section 3.2). 
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Then AC (@+ Ex) 
due to incidence, 
where Ex= (A,e/ Aja) 
and A,=Axa. 


In practice, however, the downwash distribution is not usually known with any 
accuracy even for the rigid wing, and since it will change slightly with speed due to 
wing distortion it will not be possible to carry out a single aeroelastic equilibrium 
calculation for the tailplane applying to all speeds. In view of these difficulties it 
seems acceptable to use a mean downwash angle, the mean being weighted in favour 
of the regions of greatest tailplane lift. 


Since the downwash is produced by wing lift it is logical to use an expression 
for the mean downwash angle = having the same form as the expression for C,y. 


Then € & + + An En= Eat Eyt (10) 


where the lift distribution corresponding to A°z produces a mean downwash angle 


z,, the lift distribution corresponding to C;., produces a mean downwash angle <,. 
and the lift distribution corresponding to nA, produces a mean downwash angle ¢,. 


The value of = will approximate to the value of €, defined, and 


> = Ear + Eon + Em, 
so that the use of € does not necessarily imply an approximation. 


Since the lift distributions A°z, and so on, may be obtained by using the super- 
position method, it is possible to obtain the downwash distributions corresponding 
to these lift distributions by superposition of the downwash distributions correspond- 
ing to the arbitrary lift distributions employed. This can be a somewhat involved 
process and in practice it may be acceptable to use estimated values of €., and so on, 
with semi-empirical corrections for the effects of wing distortion. 


The expression given for € is more complicated than the usual expressions for 
downwash angle, and the effect of using it in stability theory is considered in 
Section 4. 


3. Fuselage, Control Circuit and Tail Distortion 


Distortion of the fuselage, tail and control circuits modifies the tail lift and 
pitching moment contributions, and also the elevator hinge moment. This has been 
considered in Ref. 2, where compressibility effects have not been explicitly 
considered. 


Lyon and Ripley’ write 
A + A.n+A,;B8 


B, B, 2, Bn B,8 
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nr, =tail setting angle at zero wind speed (“built-in” tail setting 
angle) 


y, 8 are control angles equivalent to movements of the _pilot’s 
controls (7), 8, of Ref. 2) as is ny, if a “variable incidence” 
tailplane is used. 

Alternatively 
Cyp= A + 
Cy= + Bin + 
where 
ny= true (root) tail-setting angle 
A, r= 


= value of A, when fuselage is rigid*, 


and likewise 02+. 


3.1. Effect of fuselage inertia loading 


The expressions already given in Section 3 may be used to include the effects of 
all distortion under purely aerodynamic loading, the values of A,, A,, etc., being 
modified to take accounts of these effects. If the effect of bending of the fuselage 
under its own weight is to be included, a slight modification of the expressions for 
Cir and Cy is necessary. 


Suppose application of a (total) normal acceleration of ng produces a net change 
of nda,/0n in the value of 


0a 


where @°= true root tailplane incidence for n=0. 


This may be written 


Cir? + NA (12) 
Oay 
and similarly Ca= Cy°+nB,r 
(13) 
= B,+ + Bn +B,B. 


*This applies to tailplanes of “ fixed” incidence—with a “variable incidence” tailplane 4 ,, 
is the value of A, when both fuselage and incidence control circuit are rigid. 
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If application of ng produces a net change of tail-setting angle Ay,, we have 


A 
F; 
| where Ann= change of nr due to ng at zero wind speed 
q= 


load on tail 


F,= fuselage bending stiffness = change in yy 


Ann; 
Then Anr= 
F; 
on n on 1+A Sy 
iT F, 
4) 
qt+A, 
(Sr) 0: - =- *) (15) 
én] 


on 14 An) 


3.2. Calculation of A:, Az, etc., including distortion effects 


The analysis of Ref. 2 is adopted here but with certain modifications and with 
the inclusion of compressibility effects. 


FUSELAGE DISTORTION (AERODYNAMIC LOADING) 


Formulae are deduced in Ref. 2 for the ratios A,/a,, A./a,, allowing for 
fuselage distortion, and these expressions may be applied in general if a,, a., etc., 
are interpreted as the values of A,, A., etc., when compressibility and all distortion 
effects except fuselage distortion are included. 


TAILPLANE-ELEVATOR DISTORTION 


As pointed out in Ref. 2, tailplane distortion and elevator distortion cannot be 
considered separately, since the two are completely interdependent. For this reason 
it is not strictly possible to use the superposition method for the estimation of 
distortion effects on the tailplane-elevator combination. In Ref. 3 a method of 
treating the “ flap-deflection ” case is given which uses the superposition method, but 
this is based on the assumption that the local flap angle is not changed by distortion 
of the main surface or flap, so that the form of the spanwise load distribution due 
to the flap deflection remains the same at all speeds. While this assumption may not 
lead to serious errors when dealing with the wing-aileron combination (although it 
seems dubious even in this case), it seems likely that it would lead to appreciable 
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errors if used for a tailplane and elevator. For the treatment of the tailplane and 
elevator, therefore, the method of “semi-rigid representation” as given in Ref. 2 
would seem to be more suitable. If the effects of tab distortion and elevator skin 
distortion can be neglected (see the next paragraph), the effects of compressibility 
may be introduced (by means of linearised theory corrections) to the “strip” 
derivatives used. The corrections applied should be those appropriate to the three- 
dimensional tail surface, so that the overall lift and hinge moment coefficients vary 
in the correct way although changes in the form of the lift distribution due to com- 
pressibility are ignored, as with the superposition method. It is not possible to 
carry out the calculations in terms of a “ compressibility-distortion ” parameter such 
as qA, (used in the superposition method), since the aerodynamic coefficients a,, a., 
a,, b,, b,, b, are not all modified by compressibility in the same way. It is, therefore, 
strictly necessary to carry out calculations for each combination of g and M within 
the required range. 


ELEVATOR TAB AND TAB CIRCUIT DISTORTION 


In Ref. 2 it is shown that, if C, and C, are small (as they usually are), A,, A,, 
B, and B, are almost unaffected by tab distortion and A,, B, and C, are modified 
by the factor K /(K —C,’) where C,’ is the value of C, when tab distortion is ignored 
and K is proportional to (1/q) x tab stiffness. On this basis the elevator angle to 
trim, and hence the stick fixed static margin K,, is not appreciably changed, while 
the stick free static margin K,’ is slightly changed. Similarly the stick fixed 
manoeuvre margin H,, is almost unaffected, while the stick free manoeuvre margin 
H,,,’ is slightly changed. If the tab circuit can be designed so that K is large compared 
with C,’ for all g, tab distortion effects are then obviously negligible. If this cannot 
be done, tab distortion effects are easily included by applying the factor given to 
A,, B,, C;, on the reasonable assumption that the secondary effects of tab distortion 
on elevator and tailplane distortion are negligible. 


If power operated controls are fitted this type of distortion does not. of 
course, arise. 


ELEVATOR SKIN DISTORTION 


The effect of elevator skin distortion is also considered in Ref. 2, where it is 
shown that distortion of the elevator skin, caused by pressure differences between 
the inside and outside of the control surface, can modify A,, A, and.A, slightly and 
B,, B, and B, to quite a large extent. The treatment is approximate, since skin 
distortion of a fabric or metal covered surface is not proportional to load. 


The curves in Ref. 2 of panel deflection against load show that, for an 
unstiffened metal skin, the rate of change of deflection with load is small once the 
small initial deflection required to take any load has been exceeded. If the skin is 
applied with initial tension (as is normal practice) this small initial deflection can be 
minimised. In view of this it seems reasonable to ignore the effects of elevator skin 
distortion for design purposes on the grounds that it should be possible to design the 
elevator so as to keep the skin distortion and its effects small enough to be negligible. 
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The main effect of any skin distortion is to decrease the stick free static and 
manoeuvre margins, so that if powered controls are used the effects are even less 
serious than with manually operated controls. 


ELEVATOR CONTROL CIRCUIT DISTORTION 


The treatment of elevator control circuit distortion given in Ref. 2 can be applied 
in general if a,, a,, etc., are interpreted as the values of A,, A,, etc., when com- 
pressibility and all distortion effects except control circuit distortion are included. 


4. Modifications to existing longitudinal stability theory 


When the new expressions of Sections 2 and 3 for Ciw, Cmw, and so on, are used, 
we have, in general, 


{ A, AP — — ANE nt + NA *+A.n+A,B } (17) 
where Cur=Cm contribution from fuselage 
Vig.= Syly/ Se 


l,= distance from tailplane mean aerodynamic centre to aircraft c.g. 


The other new symbols are defined in Sections 2 and 3. 


4.1. Static margin 


In Ref. 1, 
(dn 


where CrpV? = constant = W /(4S) 
i.e. in level flight n= 1=constant. 
Following Ref. 1 the expression for C,, can be modified by the introduction of 


a new tail volume coefficient V;. Thus, when n=1, 


1 


(18) 
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ae 
*S\Aa° 
and V as. = os. + (h H,.°) Ss 


Differentiating equation (18) with respect to C, and making use of the fact that 
dn/dCy=0, gives, as in Ref. 1, 


K,= — VA, (dy/ 


Similarly for the stick free case (n= 1):— 


On 7 
(19) 
where Vac. Vac. 
1+A ) 1+F 
a 
A, 
A,= A,- B, 
A, 
2 
A, 
A,= A;- B, 


This yields similarly, K,’= —VyA, (4) 
r/¢ 
as in Ref. 1. 


4.2. Manoeuvre margin 


The manoeuvre margin, H», is given by"? 


lm, 


The second term is usually fairly small and decreases with increase of altitude. so 


that the value of the manoeuvre margin is greatly influenced by the value of 


*Note that /= distance between tail and no-tail aerodynamic centres 
294 


Diffe 


Wit 


Th 


( 
= 
= 
Th 
ap 
fo 
ta 
ve 


at 


DISTORTION ON SWEPT WING AIRCRAFT 


When n is variable, 


+ + (A (Cy Ciwe A n) (h #,") A n 


—Vy {4,[ -A,n)— )—An € ine | 


Differentiating (20) with respect to C,, with stick fixed (i.e. 1 and 8 constant), gives 


nw WwW 
If and C.,= 
Ci= nCiy and (22 (22) 


With stick free, 


The tail contribution to is, approximately, 
1 


Dn, stick fixed 


21 


The wing contribution to m, is not necessarily negligible if the wing is swept. 


In the manoeuvre margin theory it is implicitly assumed that m, is not 
appreciably changed by movement of the c.g. Thus, in deriving these expressions 
for the tail contribution to m,, the tail moment arm used is the distance between the 
tail and “ no-tail ” aerodynamic centres. This is not strictly correct since the angular 
velocity g is about the c.g. 
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When the tailplane is swept, its centre section and tips may be fairly large 
distances forward and aft of the mean tail aerodynamic centre, so that more 
accurately 


Mg tai = — Marts 


where | m,,=contribution to m, due to rotation of tailplane about its mean 
aerodynamic centre. 


In general, however, my: will be small. 


4.3. Dynamic stability—quasi-static theory 


Dynamic stability and its relation to the static and manoeuvre margins is 
discussed in Ref. 1, where it is shown that the slow divergence associated with 
K, <0 is less serious than the rapid divergence or rapid unstable oscillation 
corresponding to H,,<0. This means that the value of (0C,,/0C,)y is of greater 
importance than the value of dC,,/dC,. In this section the effects of distortion on 
the quasi-static longitudinal stability derivatives are discussed. 


If distortion effects are included on a quasi-static basis, the usual small-displace- 
ment equations of motion can be used for the flexible aircraft, but the derivatives 
have to be suitably modified. 


Ignoring the thrust contributions", 


mH), 


The largest contributions to m, and Maw, are normally from the tail, and* 


A° § A, -§ 
at (tail) ~ — A, MMq (tail)  — (25) 
These expressions are approximate only and correspond to those given in Ref. |. 
The wing contribution to maw,a is normally very small at subsonic sub-critical 
speeds, but the wing contribution to m, may be appreciable for a swept wing. An 
estimate of m, (wing) may be obtained by the method described in Appendix III. 


The derivatives x, and z, occur in the stability equations divided by », and are 
then normally small enough to be neglected. 


*For the stick free case, A, is replaced by A ~ 
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If the distortions of the aircraft components are small, it seems reasonable to 
ignore distortion effects on the drag derivatives, so that Cp, (0Cp/@a)y and 
(@C,/@M). may be estimated for the rigid aircraft. It remains to determine 
(0C,/ and (0C,,/0M), in terms of A°, A,, ete. 


Differentiating the expression for C,, given in equation (16), 


on 
Aw dn * As 


Therefore 


- Cw Cro Arr on 


(27) 


With stick fixed ($") =0, so that 


A°+ A, (1— A® és) 


A, Sr 
§ Gu 


(28) 


(= ), stick fixed = Dae 


With stick free, using equation (13), 


giving 
S: 


A°+ A, &) 
oC, 29 
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Similarly 


2%) 8% note | 
+A, |= a( A’ am ~ om 
Ais +23 om + om (3 + 
0p 
assuming (35) =(), 


and initially n=1. 


With stick fixed ( ou) 0, so that finally 


ick fixed = 


(31) 


With stick free, it can be shown that, using equation (13), 


1 0B, 1 faB, 
aM "— B, \am + am 2+ Ants) + 


(32) 


and hence finally 


sl= 


0) 
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Sy { 0A, 


én OM am"* 
stick free= — _ 
(33) 
0A, _ 9A, A, 0B, _,0A, A, 
where 0M’ 0M B, 0M ~ 0M 
TA, 2A, A, 0B, 24, A, _ 34, A,B. 24, 
Similarly 
0a M 
[ 
and 
(==) stick free= (==) +(h—H,°) A, (1— + 
da 0a 0 8. 1 
1 (a, { E 
Also 
stick fixed = aM ). + aM +a 4h + 
Oaq 0 (0a,y 
1 [ 
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2 


stick free = ( aM). {a H,°) aM A am + 


oH,” 


0A, 
+(h- He Ana Antnt no) t 


= dA°~ 4, 0% \ 
+A, 4" - am + 
I 
+ on +4" 2) aM "* B, aM 
L(G), 
“= A,(h—H,")+ Veg. [ Aye on GBD 


The value of f (initial value) to be used in these expressions will be given by 
equation (19) with C,,=0. The initial value of » may then be found from (18), again 
with C,,=0. Equation (16) then gives 2. It is assumed here that V.... is constant, 
ie. that ly is constant. This is not strictly true, since in general the mean aero- 
dynamic centre of the tailplane will move, but the error involved is probably small. 


4.4. Application of the theory 


STATIC MARGINS 
Perhaps the simplest way of determining the static margins is to use the 


formulae 
K.=-VsA, (Fr), 


Cy 0 


The values of 7 and £ to trim for a given value of Cy ~C,, are given by equations 
(18) and (19) with C,,=0. These values can then be plotted against Cp, (or M) and 
values of (d1/ and (d8/dC we =o are then obtained by graphical or 
c,,=0 

numerical differentiation. This procedure is suggested in Ref. 2, and must involve 
some loss of accuracy although, since the trim curves can be obtained by calculation, 
this loss can be minimised by the use of a sufficiently large number of ordinates. 
When the trim curves are plotted against C, the high-speed end of the range becomes 
“compressed” and the low-speed end elongated, so that it is probably better to 
plot control angles against Mach number. Then 


a =- since C,M*=constant. 
R R 
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It would be possible to find the value of dC,,/dC, from the relation 


M 


This method, however, suffers from the disadvantage that values of 0A,/0M, 
0A,/0M, etc., must be obtained by graphical differentiation (in general) so that the 
loss of accuracy is likely to be greater than that involved in finding dy/dCx. 


In general +C,=Cp Cos ye 
Ty 
Cs = Cr sin Ye 
where Tx=component of thrust force normal to flight direction 


T;=component of thrust force in flight direction 
Ye= angle between flight path and horizontal. 


Thus in level flight Cp=C, for Ty=0. This result is not, as is suggested in 
Ref. 6, dependent upon the drag being small compared with the lift. 


MANOEUVRE MARGINS 

The manoeuvre margins may be obtained by direct substitution of values for 
A’, A,, etc., in the equations of Section 4.2. It is thus a simpler matter to obtain 
the manoeuvre margins than to obtain the static margins, because of the absence of 
derivatives with respect to forward speed in the expressions for the manoeuvre 
margins. 


5. Miscellaneous Refinements 


5.1. Inclusion of effects of change of density with altitude 
on Stability 


The effect of density variation with altitude on longitudinal stability is con- 
sidered by Neumark in Reef. 18. Only the level flight condition is there considered, 


so that C, ~ Cp. 
STATIC AND MANOEUVRE MARGINS 


Neumark gives 


Cu) M 


(14 
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where N= 


and R= gas constant for air, rate. 


The static margin (stick fixed or free) as defined by Gates and Lyon") is K, (for 
C.=Crk): Kus is the “ generalised static margin” (stick fixed or free) which is pro- 
portional to E,, the last term in the stability quintic, when density variation with 
height is included. However, the quantity K,s is no longer a measure of the stick 
movement (or force) to change speed. 


From these formulae, given K, and (0C,,/0C,)™ (which is found when calculat- 
ing the manoeuvre margin) K,s can be evaluated. 


There are no effects of density variations with altitude on the manoeuvre 
margins. 


EFFECT ON STABILITY DERIVATIVES 

Additional stability derivatives xs, zs are introduced which are derived in full 
in Ref. 18. These are, however, functions of Cp, 0C)/@M, C,, 0C,/0M and of 
quantities unaffected by distortion, so that they are readily determined. 


5.2. Inclusion of effects of changes in the form of the wing lift 
distribution and movement of aerodynamic centre due to 
compressibility 


If distortion effects are to be included, it is very difficult to allow correctly for 
variations in the form of the wing lift distribution or in the aerodynamic centre 
position of the rigid wing. If they are functions of Mach number it is strictly 
necessary to perform a complete set of calculations for each combination of M 
and q. To avoid the heavy labour of such a procedure, it is suggested that if the 
shift of aerodynamic centre is not very large then the aeroelastic equilibrium calcula- 
tions can be made first ignoring the shift of aerodynamic centre. Then allowance 
for the shift can be made without correcting for the secondary distortion effects 
introduced. If the movement of the aerodynamic centre is large, however, it might 
be advisable to perform calculations with a range of representative positions, and 
then to interpolate. 


If the tailplane aerodynamic centre shift is also appreciable it may be necessary 
to modify the tail arm. 


6. A Simple Example 


As an example the hypothetical aircraft illustrated in Fig. 1 was considered. It 
has a wing of 45° sweep, an aspect ratio of 3-81, and an “all-moving” tailplane of 
similar plan form to the wing. For simplicity the structural characteristics of the 
wing and tailplane were assumed to be similar, so that calculations of lift curve slope. 
and so on, made for the flexible wing apply also to the tailplane. Compressibility 
effects were allowed for by applying linearised theory corrections to the wing 
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Fig. 1. 
Plan form of hypothetical aircraft of Section 6. 


and tail lift curve slopes, and also to an assumed value of C,,,, introduced into the 
calculations as an extra “fuselage” pitching moment contribution—the distorted 
wing has zero twist at all sections and is of symmetrical section, giving zero contri- 
bution to C,,,. No other fuselage or thrust contributions to pitching moment or 
lift were included. Tail lift was included in the total C,, but wing weight effects 
were neglected. 


Distortion of the wing, fuselage and tailplane was taken into account, but not 
control circuit distortion, and there is no tab. This is consistent with a system of 
completely rigid power operated controls, and hence the “stick fixed” case only 
was considered. 


6.1. Calculation of lift and pitching moment coefficients for the 
flexible wing and tailplane 
The first problem was to find values of A=(@C,w/02)y=(@Ciow/O%)y and 
AH, for a suitable range of values of g (=4pV*) and Mach number. The wing is an 
example of the case considered in Section 2.2.1, and hence the superposition method 
was used. 


The method used was exactly as described in Ref. 3, as follows : — 


(i) The method of Kiichemann"? was used to give (for incompressible flow) the 
lift distributions corresponding to incidence distributions e=»,, z=n,?, 
z=a constant=1 radian on the rigid wing. (),=y/{b/2}). These lift distributions 
Were integrated graphically to give values of Ci,/Ar. Ar, 
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where 


Cirz=rigid wing lift coefficient per radian of incidence 

C.,=ngid wing C;, corresponding to 2=»,, with «=1 radian at the tip 
C,,= rigid wing C,, corresponding to z=»,”, with radian at the tip 
C,;=rigid wing C;, corresponding to 2=»,*, with radian at the tip 
Ax = lift curve slope of rigid wing. 


It was found that Ci, Cr, /Axp=0°264, Cr; /Ag=0°184. 


(ii) In stage (i) the locus of aerodynamic centres was obtained, and hence the 
pitching moment distributions corresponding to «=1,, 2=7,”, and so on, were 
plotted and integrated to give the position of the mean aerodynamic centre of the 
rigid wing, and the pitching moment coefficients about that point for the various 
incidence distributions considered. 

Can _ Cus 


Then =0, —0-0775, —0-0730, 


Cas 
Ax 


where Cinz, Cm, etc., are the pitching moment coefficients about the rigid wing mean 
aerodynamic centre corresponding to Cr. C;,, ete. 


= —0-0675, 


(iii) A relation between A, and Mach number was obtained using the method 
of Collingbourne"®?. 


(iv) The lift distributions for the rigid wing were integrated to give shear 
force and bending moment distributions, it being assumed initially that the wing 
had a straight flexural axis lying along the 0:45 chord line. Torque distributions 
about this flexural axis were also obtained. 


(v) Twist and slope distributions for the four cases were obtained using 
assumed stiffness distributions, and these distributions were then modified at the 
root in an attempt to introduce corrections corresponding roughly to the root 
constraint effects on a swept wing. 


(vi) The elastic incidence changes of (v) were matched as described in Ref. 3 
and in Section 2.1.1 to give the superposition coefficients As, Bs, Cs for a range of 
values of gAr. 


(vii) Using (iii) and (vi), a graph of A and AH, against Mach number was 
produced for the condition qg/g= M7’, where g is the value of g corresponding to the 
maximum allowable equivalent air speed (Fig. 2). 


It was assumed that compressibility effects modified the two- and _ three- 
dimensional lift curve slopes without appreciably modifying the form of the lift 
distribution, or the position of the mean aerodynamic centre, of the rigid wing. 
The relation referred to in (iii) connects two- and three-dimensional lift curve 
slopes, so that when compressibility corrections are applied to the aspect ratio. 
sweep angle and two-dimensional slope, the required relation between M and A, 
is obtained. 
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Variation of A, A,, h—H, with Mach number with and without distortion. 


In (iv) and (v) it was first assumed that, with a straight flexural axis at 0°45c, 
the ratio (bending stiffness)/(torsional stiffness) was constant along the span, and 
that each stiffness varied as the cube of the local chord. On the basis of informa- 
tion given in Royal Aircraft Establishment Structures Reports 9 and 58 it was 
decided that a representative value of this ratio was 4:0, and that a representative 
root torsional stiffness was given by 
; Sl?g 
twist per unit torque at root= ee 

where /;=length of flexural axis. 


It is shown in Ref. 19 that with a swept wing of moderate aspect ratio and 
conventional construction, the concept of an effective root may be used. The wing 
may be considered to behave like an unswept wing outboard of this effective root, 
but inboard of this the root restraint effects are predominant. The information 
given in Ref. 19 suggests that for a wing of 45° sweep the effective root might be 
about 0-2 semi-span out from the root. On this basis the twist and slope distribu- 
tions were modified as shown in Fig. 11. Since this modified the overall values of 
twist and slope considerably, the root stiffness used was decreased from the value 
previously quoted, giving 
Sl?g 
“10° 

Methods of calculating the rigid wing lift distributions and the elastic 
distortions of swept wings are discussed in Appendices I and II. 


(GJ)p = 


6.2. Introduction of fuselage distortion 
Following the treatment in Ref. 2, 


A A,r 


where J=a constant, inversely proportional to the fuselage bending stiffness. 
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Using the A.P.970* fuselage bending stiffness criterion, with M=0-8 and K =0-12, 


0:05 


I= 


where gq corresponds to Vp, the “maximum allowable diving speed.” Since the 
wing and tailplane are similar, A,,.=A. 


6.3. Trim curves 


Trim curves for the example aircraft are shown in Figs. 3 and 4. Two sets 
were plotted, one corresponding to C,,,=0 (zero fuselage contribution and zero 
wing contribution about corrected aerodynamic centre) and the other to 
= —0-015/ /(1—M?) (taken as (Cino)tuseiages (Cmo)wing Temaining zero). 


It was assumed for simplicity that =.=0-1=constant, and the relation 


C,=0-0420 (4 


was used, corresponding to W/S=50 Ib./ft.? and g=4p,V?, where V = 1,000 ft. / sec. 

Curves of 7, (tail-setting angle equivalent to movement of pilot’s control) 
against C;, were produced for the rigid aircraft and for the following cases of 
distortion : — 


(i) Wing distortion only 
(ii) Fuselage distortion only 
(iii) Tail distortion only 
(iv) Wing, fuselage and tail distortion. 


6.3.1. TRIM CURVES WITH C,,,=0 

The curves are shown in Fig. 3, and it is clear that the distortion and com- 
pressibility effects introduced have had little effect on the slope of the trim curves, 
except at the highest speeds. These effects are considered in turn. 


The effect of compressibility in the absence of distortion is to displace the trim 
curve a small amount which is nearly constant for all values of C,. The curve, 
which is linear and passes through the origin when no compressibility effects are 
included, remains very nearly linear down to C,~0-1. Below this value of C.. 
the slope of the curve becomes slightly more positive, i.e. a stabilising effect occurs. 
Ref. 6 predicts that for C,,,=0 the increment of elevator angle to trim due to com- 
pressibility is very nearly constant over the whole speed range, the approximation 
becoming less exact as Mach number increases. This in agreement with the 
present results. 


*Ministry of Supply Air Publication 970. Design Requirements for Aeroplanes for the Royal 
Air Force and Royal Navy. 
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The distortion of the wing alone then produces a negligible change in the 
trim curve. The reason for this appears to be that the effect of the loss of lift curve 
slope due to distortion (the wing tips bend upwards) is offset by the forward move- 
ment of the wing aerodynamic centre. Thus, for a given C,, wing distortion makes 
it necessary to fly at a slightly higher incidence, so that for a given value of 1; the 
nose-down tail contribution is increased. This effect is almost exactly cancelled 
by the extra nose-up moment resulting from the forward movement of the wing 
aerodynamic centre, so that the value of »; to trim is unchanged. 


Distortion of the fuselage alone produces a constant increment of 1; to trim 
over the whole range, thus leaving the slope unchanged. This is because, with 
Ci» =0 and a fixed c.g. position, the tailplane load is constant over the whole 
speed range. 


Distortion of the tailplane alone also again produces a constant increment of 
tail angle to trim and the reason for this is again that the tail load is constant at 
all speeds. Since the form of the lift distribution due to any given twist distribu- 
tion of the tailplane is assumed independent of Mach number, the twist due to a 
given overall load is the same at all speeds and hence the increment of tail-setting 
angle to trim arising from twist is constant. 


When all the distortion effects are combined, the total increment of 1, to 
trim is slightly greater than the algebraic sum of the separate increments taken 
individually. This is because wing distortion makes necessary a slight increase in 
incidence for a given C,, which, for a given 7), causes an increase in the tail load. 
When the tailplane and fuselage are rigid this has little effect on »,, to trim, since it 
is largely cancelled by forward movement of the wing aerodynamic centre. When 
the tailplane and fuselage are flexible, however, the increased tail load causes 
extra tail and fuselage distortion, which in turn require a small additional 
increment of 


6.3.2. TRIM CURVE WITH C,,, = —0:-015/ /(1—M?’) 


The curves are shown in Fig. 4 and it will be seen that in this case distortion 
and compressibility effects have modified the form of the trim curves considerably. 


The effects of compressibility alone are very marked at the higher Mach: 
numbers considered. As the lift coefficient decreases and the Mach number 
increases, C,,, becomes more and more negative, i.e. the tail-setting angle to trim 
out C,,, becomes more and more negative, and at high speeds this contribution 
to nr, is large compared with that required to trim out the other wing and tail 
pitching moment contributions. At low speeds the reverse is the case, so that as 
speed increases, and the C,,, contribution becomes more dominant, the slope of the 
trim curve becomes less negative, then zero, and finally positive. 


The effect of wing distortion alone, however, then makes no further appreci- 
able difference to the trim curve, as in the case when C,,,=0. Again, the forward 
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Variation of stick fixed static margin with q/q. C,,,=0. 


movement of the wing aerodynamic centre is offset by the increased nose-down 
pitching moment contribution due to the higher tail incidence. 


Distortion of the fuselage alone makes the slope of the trim curve more 
positive (destabilising), the effect increasing with speed. The change in C.,,, 
produces a change in tail load, which becomes more and more negative (tail down) 
as the Mach number increases. As the speed increases the fuselage distortion due 
to the increasingly negative tail load produces an increasingly positive incidence of 
the tail, and therefore an increasingly negative increment of 17, is required to trim. 


Tail distortion alone has a very similar effect to that of fuselage distortion. 
The all-moving tail behaves like a wing, unlike the usual tailplane-elevator com- 
bination. Thus a positive tail load causes the tips of the tailplane to bend upwards, 
producing a positive increment of 4, to trim, and conversely for a. negative tail 
load. The changes in 7,7, to trim are therefore in the same sense as those caused 
by fuselage bending. 


It will be seen that the result of the combined distortion effects at the highest 
speeds is slightly less than the algebraic sum of the separate effects. This is 
because, as for C,,,=0, wing distortion causes a slight positive change in tail load 
which in this case at high speeds reduces slightly the magnitude of the tail load, 
which is negative, being largely determined by C,,,._ Thus the fuselage and tail 
distortions are slightly reduced and there is a small reduction in the overall 


(negative) increment of n, due to distortion when the distortion effects are 
combined. 
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6.4. Stick fixed static margin 


For an aircraft fitted with an all-moving tailplane, 


= Vid, C.=Cy). 


Values of (drir)/dC:)., 0 were accordingly obtained from the trim curves for the 
cases =O and —0-015/ /(1—M?’). 


6.4.1. 


The resulting curves of K, against q/ q for Cy,=0 are given in Fig. 5. Con- 
sider first the effects of compressibility. It will be seen that the static margin 
increases slightly with increase of Mach number. This is because (dny,/dCy),. _, 
is very nearly constant, except at the highest speeds where it becomes slightly more 
negative, while A, is increased by compressibility*. This result is in agreement 


with that predicted in Ref. 6 for the case C,,,=0. When K, is positive in incom- 
pressible flow, the restoring tail pitching moment due to a change of speed and 
corresponding change of incidence (C,V*=constant) will exceed the destabilising 
wing contribution. The difference between these contributions will be increased 
if both wing and tail lift curve slopes are increased by compressibility in roughly 
the same ratio, C,,, remaining zero. 


Considering now the effects of distortion, it is seen that, since wing distortion 
alone does not appreciably change the slope of the trim curve or the value of Vy 
or of A, for the aircraft considered, K, is almost completely unaffected by wing 
distortion. 


The effect of fuselage and tail distortion, however, is to decrease the value of 
A, (V~ still remains very nearly constant) so that, although for C,,,=0 the trim 
curve slopes are not appreciably modified by these effects, the value of K, is 
decreased progressively as qlq increases compared with the value for the rigid 
aircraft. 


In this example the loss of static margin due to fuselage distortion is greater 
than that due to tail distortion, and the combined effects produce a maximum loss 
of static margin of the order of 10-15 per cent. at a Mach number of about 0°8. 


6.4.2. Cn —0-015/ /(1— 


Curves of K,, against q/q for the case when C,,,,= —0:015/ /(1 — M?) are given 
in Fig. 6. 


Here it will be seen that compressibility effects alone produce a large loss of 
Static margin at high Mach numbers. The static margin, which at low values of M 


*V, is very nearly constant. 
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Variation of stick fixed static margin with g/g C —0°015 


mo~ 


is approximately 0-10, becomes zero at q/q-~0°5 (M~ 0-7) and rapidly increases 
negatively at higher speeds. This loss of K, with increase of M is predicted in 
Ref. 6, where it is shown that the loss of K, due to compressibility depends on the 
value of (Cino/Kn)w=o. AS shown already, the siatic margin is increased by com- 
pressibility for C,,,=0 by an amount dependent on the low-speed value of K;. 
When C,,, <0, however, the value of dC,,,/dC;, is positive and increases rapidly 
as M increases, producing an increasing loss of static margin. These two effects 
are in opposition, but the latter is dominant in this example. The loss of static 
margin will depend on the relative size of the two effects, i.e. on the value 
of / Kn) 


Considering now distortion effects, it will be seen that wing distortion alone 
causes no appreciable change in K, as in the case C,,,=0, and the reason is the 
same as in that case. Fuselage and tail distortion, however, cause a reduction in A, 
and also make the value of (dyz,/dC,), _, more positive (see Fig. 4) so that they 
produce a large loss of K,. This reduction in K,, due to fuselage and tail distortion 
is actually greater than with C,,,=0, but since here the reduction due to compressi- 
bility is very great, the distortion effects appear less important. 
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Variation of stick fixed manoeuvre margin with q/q. 


6.5. Stick fixed manoeuvre margin 
For this simple example the formula of Section 4.2 (equation (21)) becomes 


ms 4. Maw l 


The value of », was taken as 50. The wing contribution to my (Maw) Was 
estimated as described in Appendix III. 


The resulting values of H,, as functions of g/g are shown in Fig. 7. ° 


It will be seen that compressibility alone produces a loss of H,,, which steadily 
increases as M increases, and this is in agreement with the results of the analysis 
of Ref. 6. 


It can be seen how this arises from the formula 


The effect of compressibility is to increase A,, A and m,. The net result is a 
reduction in the second term, which is large, and an increase in the third term, 
which is small. 
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Comparison of values of stick fixed manoeuvre margin obtained using approximate and 
correct tail volume coefficients. 


As with the static margin, it is found that wing distortion alone has no effect 
on manoeuvre margin for the example considered. Wing distortion reduces the 
values of H, and A, so that both the first and second terms in the expression for 
H,, are increased, but the two changes almost exactly cancel each other. 


Both the fuselage and tail distortion effects, on the other hand, leave H, 
unchanged but reduce A, and m,, so that there is a resulting decrease of H,, 
which, in this example, is greater for fuselage distortion than for tail distortion. 


The combined effect of wing, fuselage and tail distortion is very nearly the 
algebraic sum of the effects taken separately. This suggests that it might be 
possible to calculate the loss of manoeuvre margin due to each and to add the 
separate contributions, but in fact this would take longer than the single calcula- 
tion for the combined effects. 


To assess the error in H,, due to using V instead of V; (i.e. due to neglecting 
tail lift in obtaining C,) values of H,, were calculated using the formula just quoted, 
but substituting (i) V..,. and (ii) V,.. for Vy. The results are shown in Fig. 8, where 
the curve corresponding to the true value.of Vy is also shown for comparison. It is 
interesting to note that the error due to using V,... instead of Vy is roughly twice 
that due to using V., and that the error due to using V.g is roughly equal to the 
magnitude of the value of m,/, (in this example). 
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6.6. Quasi-static stability derivatives 


For an aircraft with all-moving tail and no tab, the equations of Section 4.3, 
with stick fixed and taking €. as constant, simplify to 


aC, Sy 

(32), =At+ $A, 

_ 0A aA 


(Sa) ~A(h—H,)—Vee { A, = 


OCmo OH, 0A 0A 
aM aA (h- H,)- Aa —Veg (2+ A, 


Values of z, nr,, = were calculated when determining the static margin. 


The longitudinal stability derivatives involve the partial derivatives just given 
(see Section 4.3). If it is assumed, as suggested in Section 4.3, that the drag con- 
tributions are unaffected by distortion, then the effects of distortion on the stability 
derivatives may be demonstrated by evaluating (0C;,/0a)y, (OC, /0M)., (OCm/0)m, 
(0C,,/0M). with and without distortion effects. This was done for the aircraft 
considered, using the values of A, A,, H,, etc., previously calculated, and 
these partial derivatives are plotted against g/g (=M? in this case) in Figs. 9 and 
10. Only the case C,,,=0 is considered, and curves are given for the rigid aircraft 
and for the aircraft with all distortion effects included (wing, fuselage and 
tail distortion). 


THE DERIVATIVE (0C;,/02)y 
As might be expected, the curves describing the variation of (8C,/ 02)y with 


g/G q are very similar to those of A against q/ q: The tail contribution is of the 
order of 10-15 per cent. of the wing contribution, so that it is not negligible. 


THE DERIVATIVE (0C,,/02)y 

The curves for (@C,,/02)y were obtained from the values of (@C,,/0Cz)y 
calculated when determining the manoeuvre margin and the values of (0C,/02)x 
already derived. 


It will be seen that, for the rigid aircraft, (OC,,/02)y is very nearly constant 
over the whole range. Since (0C,,/0C,)y is decreased by increase of M (see Section 
6.5) and (0C,,/02)y is increased at the same time, the resultant change of (0C,,/02) 
is small. Distortion effects reduce both (0C,,/0C,)y and (0C,,/0<)y as speed is 
increased, so that (@C,,/@2)y decreases steadily with increase of M for the flexible 
aircraft, 
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THE DERIVATIVE (0C;,/0M). 

Since (0C;,/0M). is a “constant incidence” derivative, the expression for it 
contains terms of the type a0A/0M. The magnitude of (0C,/0M). therefore 
depends on the magnitudes of 2, a, and 0A/0M, 0A,/0M. For small values of 
q/q, 2A/0M is small but « changes rapidly, while for high values of g/g, 2A/@M 
changes rapidly but « then changes slowly. The net result is (@C,,/@M). has a 
maximum value for moderate values of g/q. 


Distortion increases 2 and reduces 0A/0M; the latter effect is dominant and 
there is consequently a reduction in (0C;,/0M). due to distortion which is nearly 
constant over the whole range. 


THE DERIVATIVE (0C,,/0M). 


Since 0C,,,/@M=C,,,=0, and for the rigid aircraft (k—H,) is constant, i.e. 
0H,/0M=0, therefore the expression for (@C,,/0M). quoted is reduced to two 
terms : — 

0A ~ 0A, + 0A } 
The second term is small compared with the first (i.e. the wing contribution is large 
compared with that from the tail) so that 


H,). 


Since (h--H,) is constant, the form of the curve of (0C,,/0M). for the rigid 
aircraft is very similar to that of (@C,/0M)., and thus 


OM /, (h- H,) 


When distortion effects are introduced, 0H,/@M is no longer zero, and in fact 
the term (—A20H,/0M) is dominant, being large and positive. At low speeds, 
where a increases rapidly as q/q is reduced, this term causes a large increase in 
(0C,,/0M).. At higher speeds the increase is somewhat less although still 
considerable. 


EFFECTS OF DISTORTION ON STABILITY DERIVATIVES 


For the aircraft considered, with C,,,=0, it is seen that distortion reduces 
and (@C,/@M). so that Zw, Mwy and Z, are decreased. The 
effects on Zz, and my are appreciable even at moderate speeds, but the effect on Zz. 
will be serious only at high M and low C,. Distortion increases (0C,,/0M), so that 
m, is increased, and this effect is most marked at the lowest speeds. Tail contri- 
butions to m,, Maw,at Will be decreased by distortion since both A and A, are 
reduced. On the assumption that drag derivatives are unchanged by distortion, 
X, and xy will not be affected. 
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It must be emphasised that the foregoing results apply only to the aircraft 
considered and for the condition C,,,=0, and caution is needed in attempting to 
generalise from them. 


7. Concluding Remarks 


There is no intrinsic difficulty in including distortion effects on longitudinal 
stability if frequency effects are neglected. The amount of calculation involved 
may be considerable, however, especially if wing twist, camber and weight effects 
have to be included. If, therefore, there is good reason to believe that distortion 
effects are small (as they may be on a fighter-type aircraft with high stiffnesses) a 
simple crude assessment may be adequate. On a large bomber or transport aircraft, 
for which load factors and hence stiffnesses are likely to be lower, and which are 
likely to have higher aspect ratio wings, it may be essential to make a detailed 
analysis of distortion effects on the lines discussed. This is particularly true of the 
“nodded engine” layouts used in conjunction with thin swept wings. 


For the aircraft considered in Section 6 the results show that if C,,,=0 neither 
compressibility nor distortion effects modify the static margin to any serious extent. 
When C,,, is negative, however, both effects are heavily destabilising as far as the 
static margin is concerned. It therefore seems advisable to give a high-speed 
aircraft a layout which is as nearly symmetrical as possible about the plane 
of the wings. 
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APPENDIX I 
CALCULATION OF AERODYNAMIC LOADING ON RIGID SWEPT WINGS 
Several methods of calculating the lift distribution on swept wings have been 
produced. Among these are those due to Multhopp,®? Garner,“ Falkner, Weis- 
singer,"?) De Young,“?) Kiichemann ‘” and Diederich."*) Stanton-Jones@® has derived 
empirical formule based on the results of Weissinger which enable lift distribution to 
be calculated very rapidly for trapezoidal untwisted wings. None of these methods 
includes thickness, viscosity or compressibility effects. 


Of the methods quoted, those due to Falkner, Multhopp and Garner might be 
described as “lifting surface” theories, by means of which estimates of span and 
chordwise lift distributions, lift curve slope and aerodynamic centre locus position may 
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be obtained. All three methods require a considerable amount of computation if the 
full advantages of the “ lifting surface ” method are to be obtained. 


The other methods quoted are mainly modifications of the Prandtl “ lifting line ” 
theory used for unswept wings and suffer from the disadvantage that reliable estimates 
of the aerodynamic centre locus cannot be obtained. However, they require less 
computation than the true lifting surface methods. The method of Kiichemann is 
different from the others in that it uses assumed chordwise lift distributions for the tip, 
approximately mid-span, and root sections, with distributions at the intermediate sections 
based on an empirical interpolating relation. A modified version of the flat plate 
“loading law” is used. It is thus possible to obtain an estimate of the aerodynamic 
centre locus, although the estimate cannot be more accurate than the initial assumptions 
involved. The method is, however, simple to apply and provides a reasonably accurate 
estimate of the spanwise lift distribution, and it was used in the example of Section 6 
in this paper. 

APPENDIX II 
CALCULATION OF DISTORTION OF SWEPT WING UNDER GIVEN AERODYNAMIC LOADING 


The mode of distortion of a swept wing under load is different from that of an 
unswept wing, because of the effects of the oblique restraint at the root. Whereas an 
unswept wing can be treated as a simple cantilever beam in bending and as a tube in 
torsion, if the concept of a straight flexural axis is accepted, a swept wing presents a 
more difficult problem and the flexural axis concept cannot strictly be used. Near the 
root the wing tends to bend about the root (streamwise) chord, due to the restraint, but 
the outer portion of the wing behaves more like a straight wing, i.e. it bends about a 
line roughly normal to the leading and trailing edges. Just outboard of the root there 
is a “transition” between these two types of distortion, so that the overall mode of 
distortion is complicated. 


To estimate deflections for aeroelastic calculations it is possible to use the simple 
beam theory for a high aspect ratio swept wing. The root restraint effects are confined 
to a small fraction of the semi-span when the aspect ratio is high. The wing is then 
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treated as a straight wing from the structural point of view, but the slopes and twists 
are resolved into the line of flight to give the incidence changes. There is some 
ambiguity about the position: of the effective root, and this method can only be 
approximate. It is, however, quite widely used and it is put forward in Refs. 3, 20 
and 21. 


Experimental evidence''®) shows that as the aspect ratio is decreased and the sweep 
angle is increased this “ effective root” method ceases to yield useful results, because 
the region in which root restraint effects are large now occupies a large fraction of the 
semi-span. The deflections of quite a large portion of the wing cannot then be estimated 
by this simple approach. It is possible to make a correction for moderate aspect ratio'° 
but for low aspect ratios (or where root stress distributions are required) it is necessary 
to use a more refined approach. The simple beam theory can be applied together with 
a “ self-equilibrating ” stress system,'*?) and this method has been used successfully in 
practice. More exact methods have been put forward (e.g. Ref. 23) but these are 
somewhat more complicated. Once the simple beam theory is discarded, the simple 
integrations previously used to obtain deflections can no longer be used, and matrix 
or some other methods must be employed in conjunction with influence coefficients. 
A review of recent methods considered in the United States is given in Ref. 24. 


Finally, it must be noted that with a very low aspect ratio (delta) wing, there may 
be changes of camber which produce effects of the same order of magnitude as the 
twist effects. 

APPENDIX III 


ESTIMATION OF WING CONTRIBUTION TO ™, 


With a swept wing, the wing contribution to m, may be appreciable compared with 
the tail contribution. For the example of Section 6, m,,, was estimated as follows. 


Following Ref. 1, a change of manoeuvre margin is obtained for an increment of 
wing pitching moment AC,,,,, arising from the wing contribution to m,. This is given by 


AC, is the pitching moment increment due to rotating the wing with constant angular 
pitching velocity g. Suppose all the wing lift to be generated at the quarter-chord line, 
and that the wing is rotating about the point A, distant x, from the root quarter-chord 
point taken as the origin. 
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The incidence change at the point (x, y) on the quarter-chord line due to q is then 


— 4 
Aa=+ 


x=y tan 45°=y, 


The component —qx,/V is constant along the span. 


For the purpose in hand it is sufficient to neglect distortion effects, i.e. may is 
determined for the rigid wing. 


The constant incidence change — qx,/V provides zero pitching moment contributions 
about the rigid wing mean aerodynamic centre. Therefore only the incidence change 


need be considered. This is an incidence distribution that increases linearly with span- 
wise distance from the root. 


The calculations of (ii) in Section 6.1 for the case z=n,, give C,,,/Agz=—90-0775 
per radian incidence at the tip. Hence the pitching moment coefficient about the rigid 
wing mean aerodynamic centre due to the rotation at angular velocity q is 


AC = — 00775 


In this case, 


Also 


0:0775 Ay X 0-907 
Therefore =+ 


The corresponding tail contribution is +VA,/(2;). A, and Ap are of the same order. 
and V,, ~ 0-3, so that in this case 0°25 

Note that AC,,,, is the pitching moment about the mean aerodynamic centre and 
not the centre of gravity. This is consistent with the use of / instead of J, in the 
formula Aa;=ql/V and also with the assumption already made in Section 4,2 that 
Myy is independent of c.g. position. 
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